"Studies on the passage of heavy- 

CHARGED PARTICLES THROUGH MATTER 


A Thesis Submitted 

In Partial Fulfilment of the Requirements 
for the Degree of 

DOCTOR OF PHILOSOPHY 


By 

BRIJESH KUMAR SRIVASTAVA 


to Che 


DEPARTMENT OF CHEMISTRY 

INDIAN INSTITUTE OF TECHNOLOGY KANPUR 

NOVEMBER, 1974 



i 


STATSMBNT 

I hereby declare that the work presented in this 
thesis is the result of investigations carried out by me 
in the Department of Chemistry, Indian Institute of Tech- 
nology, Kanpur, India under the supervision of Professor 
Shankar Mukherji. 

In keeping with the general practice of reporting 
scientific observations, due acknowledgement has been made ! 
wherever the work described is based on the findings of 
other investigators. | 

: . r A ^ . ■! 

i 1 1 .-V ■ O' 1 1 
Brijesh Kimar Srivastava 







11 


DSFARTMSNT OF CHEMISTRY 
INDIAN INSTITUTE OP TECHNOLOGY, KAU3PUR 

CERTIFICATE I 


This is to certify that Mr. Brijesh Kumar Srivastava 
has satisfactorily completed all the courses required for 


the Ph.D. degree programme. These courses include: 


Chm. 50 1 
Chm 521 
Chm 523 
Chm 524 
Chm 541 
Chm 542 
Chm 543 
Chm 600 
Chm 622 
Chm 800 
Chm. 80 1 
Chm 900 


Advanced Organic Chemistry I 
Chemical Binding 
Chemic al Hie rmodyn amic s 
Modern Physical Methods 
Advanced Inorganic Cliemistry I 
Advanced Inorganic Chemistry II 
Nuclear Chemistry 
Mathematics for Chiemists 
Chemical Kinetics 
General Seminars 
Graduate Seminars 
Post-Graduate Research 


Mr. Brijesh Kumar Srivastava was admitted to the 


candidacy of the Ph.D. degree in December 1971 after he 
successfully corrpleted the written and oral qualifying 
examinations. 

(A. ChaKravorty) f' 

Head/ 

■ Depaartment of Chemistry 



(P.-R. SINGH) 
Convener 

Departmental Post- 
Graduate Committee 



CBRTIFICATB II 


Certified that the work presented in this thesis 
entitled: "STUDIES ON THE PASSAGE OF HEAVY CHARGED 
PARTICLES THROUGH MATTER" has been carried out by Mr. 
Brijesh Kumar Srivastava mder my supervision and the 
same has not been submitted elsewhere for a degree. 


Shankar Mukherj i 
Thesis Supervisor 





iv 


AGKNOWLEDQBMBNTS 


I am highly indebted to Professor Shankar Mukherj i for 
his valuable guidance, inspiration, and encouragement through- 
out the course of my research work. I should like to be espe- 
cially grateful to him for helping me to develop increasing 
awareness and insight not only in the area of research acti- 
vity, but also in- various other fields of general interest. 

I am highly thankfxil to Professor L.C. Korthclif fe, 

■Texas A & M University, Professor L. Bridwell, Murray State 
University, for kindly providing original data from their log 
books, and Dr. C.D. Moak, Oak Ridge National Laboratory for 
sending some of his old piiblications . 

I am grateful to my colleague Dr. Tej asvi Sharma for 
his patient perusal of the thesis. 

I wish to thank my colleagues and friends messrs. Kamal, 
Beni, Ashok, Naidu, Prabhakar, Manmohan, Vijay, Drs. Prakash and 
Kothari for their cooperation during the research work. 

I wish to express my thanks to the Department of Atomic 
Energy, for financial assistance in the form of a fellowship. 


Bj^-ijesh Kumar Srivastava 



PEBFAOB ■ 


V 


Ihe thesis entitled, ‘ Studies on the Passage of Heaw 
Charged Particles Through Matter', deals with the formulation 
of range-velocity and stopping power equations for fission 
fragments, partially stripped and completely stripped heavy 
ions in solid as well as in compound media and the calculation 
of the effective charge (mean equilibriiim charge) of a moving 
ion . 

Chapter I introduces ion penetration phenomena and 
describes broadly the need of the present investigation. 

Chapter II describes briefly the range and energy loss 
concepts and some experimental methods for the determination 
of ranges and energy-losses.. 

In Chapter III various theories, such as the classical 
theory of Bohr, Bethe' s quantum mechanical approach, Bloch' s 
formulation, and the semi-quantum mechanical treatment of 
Bohr, for the study of interaction of charged particles with 
matter are described. 

Chapter IV describes two universal range-velocity and 
energy-loss equations for fission fragments and partially 
stripped heavy ions in solid and compound media, obtained 
through the present investigation. It is found that in the 
case of light media (Z<45.5), the number of orbital electrons 



n(U ) having velocity less than a given velocity U is given 
S' ' s 

by^ 

n(Ug} =0.28 , 

o 

2 

where Z is the atomic number of the mediijm and V = e /h. 

o 

Chapter V deals with the development of the stopping 
power equations for completely stripped heavy ions. Ihese 
equations have been obtained by considering carefully the 
limits of integration of the logarithmic terms in Bohr’s 
energy loss equation. 

In Chapter VI, calculated values of the effective 
charges of light ions (Z^S.S) in gaseous media are presented 
and are compared with the corresponding experimental values. 
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CHAPTER I 


INTRODUCTICaJ 

The passage of energetic heavy charged particles 
through matter still continues to be of great interest 
even after sixty years of active investigations.^ This 

interest arises out of several theoretical and practical 

. . 2 
considerations .In 1913/Niels Bohr laid down the theoreti- 
cal principles for solving this problem on the basis of 
coulorribic interactions between the penetrating charged 
particle and the orbital electrons of the penetrated medium. 
Bohr' s treatment involves the notion of impact parameters 
and his great contribution lay in postulating a definite 
"cut off" point for the impact parameter and in providing 
a simple, but approximate value for the maximum impact para- 
meter. In 1948, Bohr reconsidered the whole problem and 
went to great lengths in examining the limits within which 

the classical orbital picture and the idea of impact para- 

■ 4 

meters may be assumed to be valid. Meanwhile Be the had 
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developed his quantum mechanical treatment for calculating 
the energy loss of charged particles in matter and Bohr^ 
proposed a generalized stopping power, formula which incor- 
porates the features of his earlier classical formula^ as 
well as the features of Bethe* s quantum mechanical formula.^ 
However, when one attempts to use either Bethe* s formula^ 
or Bohr's formula, one comes across many practical diffi- 
culties . 

Firstly, there are definite domains determined by the 
nature of the penetrating particle and the penetrated media 
(mass niomber and nuclear charge) and the velocity of the 
particle, within which these fomulae are valid. Apart from 
these restricted domains, there are broad regions where none 
of the above treatments are applicable. For practical purposes 
one often needs compute the energy-loss of particles which 
belong to these "no man' s land". Further, even within their 
separate restricted regions of validity, each formula cannot 
be directly used for stopping power calculation because each 
contains some parameters whose values are not obtainable 

through rigorous theoretical calculations. Ihus, inspite 

• 2 . 
of the existence of Bohr's earlier formula for sixty years 

and Bethe Vs formula'^ for forty years, most of the practically 

used stopping power formulae suffer from considerable amount 

of empiricism. ^ Secondly, rapid development has taken place 
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in the generation of mediian energy charged particle accelera- 
tors in the last decade. This development has revived 
interest in the stopping powers for the heavy ions produced 
by these accelerators. In the current literature termino- 
logy, "heavy" ions mean particles much heavier than protons. 
These heavy ions, unless they are accelerated to a very high 
energy, remain partially stripped of their orbital electrons, 
and while passing through matter gradually pick-up electrons 
and become progressively neutralized as they slow down. This 
phenomenon introduces a further complicating factor in the 
stopping power formulae mentioned earlier, since these were 
based on the assumption that the ionic charge state of the 
ion remains invariant throughout. A majority of these heavy 
ions belong to the domain in which neither Bohr’ s classical 
treatment nor Bethe's quantum- mechanical formulation is 
valid. Fission fragments belong to this class of heavy 
ions and because of their practical importance in a variety 

of fields such as fission -chemistry^ and direct energy conver- 

6 

sion in fission-electric cells,' the total penetration depths 
or ranges of fission fragments as well as their rate of energy 

loss in different materials, have been the sxfDject of nume- 

^ . 7-10 

rous investigations. 

Recently, Lindhard et al.^^ put forward a stopping 
power formula which v/as assumed to hold for heavy ions and 
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fission fragments. However/ disagreenrent exists between the 

experimental data and the corresponding theoretically calcu- 

8 

lated values. Practical range-velocity relationship remain— 

12 

ed completely empirical till Niday proposed hxs range-velo- 
city equation. This equation was obtained by modifying Bohr* s 

3 

generalized stopping power formula and was found to be fairly 
precise in the case of a large nurrber of fission- fragments 
penetrating metallic uranium. The purpose of the present 
investigation is (a) to obtain very general stopping power 
and range-velocity equations for partially stripped heavy 
ions and (b) to derive stopping power equations for completely 
stripped ions which belong to neither of the domain*^ in which 
Bohr' s or Be the' s equation is valid. 




CHAPIER II 


EXPERIMENTAL METHODS FOR DETER- 
MINING THE VELOCITY, RANGE, AND 
ENERGY-LOSS OF FISSION FRAGMENTS 
AND CHARGED HEAVY PARTICLES 

II. 1 FISSION-FRAGMENT VELOCITY 

1 O 

The tinie-of-f light technique, originated by Leachman 
consists of finding out the time required by a fragment to 
travel a fixed distance from its point of origin and thus 
its velocity can be obtained. The time of origin is provi- 
ded by the pulse from the electrons emitted from the backing 
of the fissile source travelling the 1 cm. distance from the 
fission source to the nearest scintillation detector. The 
time-of-f light of the complementary fragment through the 
343 cm drift distance determined the time of occurrence of 

the pulse from the remote detector. Utilizing this technique 

13 , 

Leachman measured the velocity distribution of fragments 

f rom the f ission of and ^^^Pu. Comparison of the 

results of the velocity distribution from the time-of- flight 

technique and ionization measurement shov7s that the latter 
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method gives lov/er values for the kinetic energy due to an 
inherent * ionization defect' . 

.14 

Stein used two evacuated drift tubes with a thin 
fissile material placed between them, at the end of which 
two scintillation detectors are placed. Ihe zero time is 
recorded from a pulse generated by a large number of electrons 
(^50- loo) stripped out of the thin backing foil, when one of 
the fragments passes through it on its way down the drift 
tube. Ihese electrons are accelerated to a high potential 
and focussed on the photomultiplier tube, the signal produced 

thus gives the beginning of the fission event with a time 

. -9 

resolution of 5 x lo sec. 

II. 2 RANGE-CONCEPT AND IMETHODS OF RANGE-DETERKENATION 

Ihe total distance travelled by a charged particle in 
an absorbing medium before it comes to rest is equal to its 
range. Asstming a laniform linear irotion, the greater 
kinetic energy, the more would be the penetration depth. A 
number of techniques have been used for range determination 
of the fission products. 

II. 2.1 Fission Products* Ranges 
(a) Thin target method: 

A target can be thin or thick as compared to the range 
of the particular recoiling product in the target. For example. 
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2 

for a range of 20/^g^cm of the particular recoil a target 
2 

of 1 Ag/cm is a ^ thin target' . If the source is very thin 
and the collimator is so narrow that the recoil products tra- 
velling only at right angles to the catcher are allowed to 
penetrate the catcher then, theoretically, all the products 
of identical mass should come to a stop at a particular depth 
of the catcher, which is equal to the range of the product in 
that catcher material. However, this does not happen in 
experiment because of the combination of two factors: (a) on 
account of the statistical nature of the collision between 
the products and the absorber material, which determines the 

CO 

energy loss of former, some products come totstop a little 
before and some a little beyond the thickness at which they 
should stop and (b) there is an initial spread in the kinetic 
energies of the products of identical mass . The ranges with 
thin targets are determined by two miethods. 

(1) Thin target-thin catcher differential method ; In 
this method, a well collimated beam of fission products, emer- 
ging from a thin layer of the fissile material is allowed to 
penetrate the absorber xvhich is either in the form 'of stack 
of thin foils or a thick foil from which thin layers are 
removed later either by chemical stripping or by fine grind- 
ing. The activity in each thin foil or each layer removed 
is measured by radiochemical techniqj.es. 
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From the plot of the activity in each layer vs the thick- 
ness of that layer one may obtain the mean or most probable 
range of a particular fission product. This method requires 
sufficient activity in each thin layer so that counting errors 
may not be too high. The possible sources of error in this 
method lies in the finite thickness of the fissile material/ 
the finite size of the collimating holes, and in the non-uni- 
formity in the thickness of the thin foils or the thin layers 
removed. 

(2) Thin target-thin catcher integral method: In this 

method a thin target and several thin catcher foils are stacked 

together and irradiated for a long time without using a colli- 

mater. The activity 'A^' present beyond a thickness * t' from 

the surface of the stack facing the target is related to the 

17 

range ' R' as given by Segre and Vjiegand^ 

A^ = C (R-t) , 

where C is a constant. Thus, a plot of A^ vs t gives a 
straight line the extrapolation of which to the ' t' axis 
gives the range R. 

(b) Thick Target Method :. 

Thin target methods become unsuitable for low fission- 
cross section materials, particularly when fission is induced 
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by accelerated charged particles or fast neutrons. The main 

disadvantage with accelerated charged particles is that the 

beam intensity decreases with increasing energy of the pairticle 

and to have a measurable activity^ the thin target-thin 

catcher method would require long periods of irradiation. 

18 

Sugarman et al. developed the thick-target integral and 
differential methods for the range measurements. The inte- 
gral method gives the range value in the target material 
while the range in catcher material is obtained by the diffe- 
rential n^thod. 


1 2 

(1) Thick target-thick catcher method ; Niday used 

this method for measuring ranges of nuclides from the thermal 

235 

neutron induced fission of U. A thick foil of the target 
material is sandwiched between two absorber foils (thick 
enough to stop all the recoiling nuclides) and irradiated 
with thermal neutrons. If it is assumed that the fission 
fragments are emitted isotropically, then the range R of the 
nuclide in the target material is given by - 


A, 


R = 




4 . t 


where A is the activity of the nuclide in one of the catcher 
C 

foils, A^ is the activity of the nuclide in the target and 
* tv is the thickness of the target foil. Chinaglia et al. 
have used modified form of this method to determine the 
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r'a,ng©s in cairbon using high resolution Ge(Li) cietector* With 
this technique it has been possible to measure the ranges of 
several fission products without doing any chemical separa-' 
tion. A homogeneous source was obtained through thermal diffu- 
sion of enriched uranium at a high temperature into pyrolytic 
carbon discs* But this technique is limited to nuclides that 
have reasonable yields and emit prominent gamma rays., Gordon, 
et were the first to develop this technique. Ihe 

accuracy of range measurements obtained with Ge ( Xi) spectro- 
metry is potentially at least as high as for the standard 
radiochemical methods, and the former can be used for volatile 
products that are difficult to handle radiochemically. 


( 2 ) fhlck target-thin catcher method ; Alexander and 
Gazdik^ used moderately thick targets sandwiched between 
stacks of absorber foils of different thicknesses. Ihe average 
range R of a product may be calculated from the equation: 



t _ -g-W. 
R 2R 


where is the activity of the recoils penetrating beyond a 
thickness ‘ t‘ of the catcher assembly, R is reciprocal of the 
average range, w is the target thickness and C is a constant 
for converting the rate of velocity loss in the target mate- 
rial to that in the catcher material and is evaluated from 
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the expression: 

U 

A w 


C_ 

4R 


11 *2.2 Ranges of Accelerated Heavy Ions 

The differential range of accelerated heavy ions is 
obtained by performing integration of the energy- loss data: 

R (E) - R (E ) = f dE 

o J dx 

Eq 


where R(E) and E.(E^) are the ranges corresponding to energies 

E and E^ respectively. Using this procedure Bridwell and 

09 79 127 

Moak'^'^ determined the ranges of Br and I ions in Be, Al, 

Ni, Au and UP^ . 

11.3 ENERGY- LOSS MEASUREFJSNT 

Experimental measurement of the stopping powers of heavy 
ions requires accurate knowledge of the energies of the ion 
before and after it has passed through a thin absorber foil of 
precisely known thickness. The ratio of the energy loss to the 
foil thickness is regarded as the stopping power at a mean 
energy defined by , , 

E = E - , 

m 1 2 

v/hefe E^ is the incident beam energy, b^E is the energy loss 
in foil and 
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II. 3.1 Fission Fragments 

Hie experimental set-up for the m.easurement of energy- 

loss of fission- fragments consists of absorbers of different 

thicknesses mounted on rings and positioned approximately on 

the mid-plane between the fissile source and the solid state 

detector in an evacuated chamber. If only one detector is 

used, then the energy-loss of the average fission fragment can 
23 

be obtained. Using two detectors on both the sides of the 

fissioning material ( thin-source) , the energy loss of the 

median light and the median heavy fragments may be obtained 

9 24 

by appropriate coincidence circuits. But none of 

these measurements give information regarding the energy loss 
of each of the observed fission fragments (variation of .stopp- 
ing-power with mass) . Hie experiment for the measurement of 

the energy loss of individual fission fragments involves the 

lO 

evaluation of three parameters : the energy of the primary 
fission fragments, the time necessary for it to traverse a 
measured distance and the energies of the corrplementry fragments. 

The fission fragments after passing through the absorber, 
impinge on the detector and the recorded pulse-beight spectrum 
yields the fragment energy as follows: 

25 ■■ 

Schmitt et al. showed that the energy E of the 
ment is related to the pulse— height and the mass A of the 
fragment by the expression: 
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S = (a + a' A) P + b + b’ A 

where a, a' ^ b and b' are detector constant and given by 

a = 24.0202/(P^-Pj^) 

, a' = 0.035 74 / (Pj^-Ph) 
b = 89.6083-a P^^ 
b' = 0.l370-a' Pj^ 

Pj^ and are the peak positions of the light and the heavy 
fission f ragmen tsv and P is the pulse height as determined 
from a fission fragment kinetic energy spectrum, 

II. 3. 2 Charged Particles 

A beam of accelerated heavy ions is passed through an 
analyzing magnet of high resolution to ensure that the initial 
energy is constant and well defined. This mono- energetic beam 
is then allowed to travel a known thickness of the absorber 
foil placed in its path. The energy of the transmitted ions 
is then measured using either an electrostatic analyzer equi- 
pped with a position sensitive detector, a magnetic spectro- 

graph,^^ a solid state detector^^ or a time-of -flight equip- 

.26 
ment , 

The measurement of the absorber foil thicknesses is done 
either by weighing each foil on a semi-microbalance and 
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measuring hhe diameter of foils with a comparator or by the 

'^-particle energy-loss method. In the latter procedure the 

thickness of each absorber foil is determined by measuring 

252 244 

the degraded energy of the Cf or Cm cc-particles. Use 
is made of the ci-particle range and stopping power data. 



CHAPTER III 


THEORETICAL ASPECTS OF 
STOPPING PROCESS 


111.1 ENERGY-LOSS PROCESS 

The mechanism by which a swiftly moving charged particle 
loses its kinetic energy involves four types of interactions. 

111. 1.1 Inelastic collision with the atomic electrons 
Inelastic collisions involve excitation and ionization 

of the atoms of the medium.. The energy imparted by the moving 
particle to the atoms of the medium causes ionization. Presxim- 
ably, both ionization and excitation occur simultaneously. The 
slowing down of the moving particle due to this phenomenon is 
also known as "electronic stopping". 

111. 1.2 Elastic collision with the nucleus 

This type of interaction / also described as nuclear 
collision, is equivalent to a two-body billiard ball like 
collision between the moving particle and the heavy nucleus 
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of the stopping mediiim. This process involves the transfer 
of kinetic energy to the stopping atom as a whole. 

111. 1.3 Inelastic collision with a nucleus 

A close non-capture encounter of the moving particle 
with the nucleus of the stopping atoms causes a deflection 
of the moving particle and frequently leads to the emission 
of a quantum of radiation appealing in the form of bremsstra- 
hlung. The intensity of the radiation is inversely propor- 
tional, to the square of the mass of the particle and is not 
of much importance in the case of heavy ions. 

111.1.4 Elastic collision with the atomic electrons 

In this type of interaction the elastic deflection 
of the moving particle in the field of the atomic electrons 
of the stopping atoms takes place. The energy transfer is 
generally less than the lowest excitation potential of the 
electrons. Such collisions are significant only in the case 
of very low energy electrons. 

Electronic collision is usually the predominant 
mechanism by which a charged particle loses its energy. 
Nuclear collision starts dominating when the ion velocity 
becomes less than the velocity of the electron in the first 
Bohr orbit. In the following section the theory of electronic 
collisions will be presented. 
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III, 2 ELECTRONIC STOPPING PROCESS 

A detailed discussion of the interaction between a 
moving charged particle and an electron of the stopping medixom 
was put forward by Bohr^ in 1913. Using classical mecha- 
nical ideas an ej<pression for the rate of energy loss was 
obtained. Later on in 1930 Bethe derived an expression for 

the energy loss of high speed light particles. Three years 

29 

after Be the' s quantum mechanical treatment^ Bloch succeeded 
in getting a stopping poxver equation which incorporates the 
classical treatment of Bohr ^ and the quantum- mechanical deri- 
vation of Beths'^ . In 1948 Bohr^ again considered the pheno- 
menon of charged particle energy-loss and derived a stopping 
power equation which incorporates both classical as well as 
quantum mechanical ideas. 


III. 2.1 Classical theory of Bohr 

The interaction between a charged particle and an 
atomic electron is described as a coulomb interaction between 
the electron and a positive point charge. Thus, iftheavy ion 
decreases its charge by picking up electrons in its passage 
through matter the rate of energy- loss will diminish. The 
energy transfer to an electron in a collision with an ion of 
velocity V is given by 



o 2 4 
2z e 


2,-2 

m p V 



. . . (3.1) 
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where m and e are the electronic mass and charge/ p is 
defined as the closest distance of approach of the moving 
ion to the electron and is termed as 'impact parameter' and 
ze is the charge of the moving ion. In the derivation of 
Eq.(3.1) it is assumed that 

(a) the electron is free and initially at rest, because 
in each collision the struck electron receives so much energy 
that its initial binding energy can be neglected. This is 
justified as long as the velocity of the ion is greater than 
that of the electron. 

(b) The ion moves almost undeflected, because the mass 
of the electron is negligible conpared to that of the heavy 
ion. 


In Eq.(3.1), AE is the energy gained by one electron. 
The energy gained by the electrons lying between iirpact para- 
meters p and p+dp and in a small thickness dx of the stopping 
medium is given by 

(3E = As X number of electrons = AE . 27fp n op dx 


dE =, 


2 4 

47rz e n 

y dx dp 


... (3.2) 


m p V 

where n is the number of electrons per unit volume. The total 
energy loss per xinit path length is given by 


dS _ 4TTz^e^n \ - 

dx ~ „2 ) p 

m V r' . 


... (3.3) 
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The limits of integration are given by the values of p 

insLx 

and p^.„. The maximum value for the impact parameter p 

min . jr i' -^max 

is given by the time of interactionT of the ion with the 
electrons, which should be shorter than the natural period 
of vibration of the electron. The upper limit for the inpact 
parameter is then 


max 


= Tv 


1.123 V 
CaI 


. . . (3.4) 


where Oo is the cyclic frequency of the electron in the atom 
of the medium. The mdnimum inpact parameter is determined by 
the maximum amount of energy that can be inparted to ah elec- 
tron. From Eq.Cs.l) 

2 

Prr,-iv. = ...(3.5) 

min m 

Using Eqs.(3.4) and (3.5) in Eq. (3.3), the total energy loss 
becomes 


dx 


2 4 

4^2 e n 
m V 


In 


1.123 m 

2 ' 
z e to 


. . . (3.6) 


the 

if all the electrons have^same characteristic frequency 6 o 

and ^.^ .is the atomic number of the medium. If each atom 

contains Z 2 electrons having individual frequencies 6 U 2 , • • 

(jj , the average energy loss becomes 
s 


47 r 2 ^e^n ^^^2 1.123 m 

m s=l 2 e ^ 60 ^ 


dE 

dx 
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Sq.(3.7) was first derived by Bohr^ using sirtple classical 
mechanical ideas. 

III. 2. 2 Criterion for the applicability of the classical 
orbital picture 

Later in 1948 Eohr^ showed that the classical theo3:y 
is valid only for those collisions in which the ion velocity 

V greatly exceeds the orbital velocity U of the atomic elec- 

s 

trons in their respective orbits. This imposes the restric- 
tion that V should exceed the velocity of the inner shell 
electrons (K -Shell) for ionization in the K -Shell. The 
second difficulty arises due to the fact that the maximum 
inpact parameter must exceed the radius of the Bohr orbit 

(p \ a ) . This gives the minimum energy transfer to the 
mox/ s 

electrons, which are less than the ionization potentials and are ^ 
therefore, inconsistent with classical ideas regarding the 
energy transfer to orbital electrons. The impact parameter of a 
particle approaching a force centre can be precisely deter- 

■“ i 

mined only v/hen the incident particle is well collimated. 

This collimation can be achieved only with a narrow 4perture 
but which would, however, cause a diffraction of thd particle. 

By considering the diffraction caused by the aperture and t^ 

3 , 

deflection of the particle, Bohr showed that 

%)} Ir .. . (3.8) 
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the ■ 

is the essential condition for the validity of^classical 

2 

orbital picture, where •%= 2 ze/(ilV) and ze is the charge of 
the moving particle of velocity V. For decreasing values of 
the orbital picture gradually loses its validity. 

III. 2. 3 Bethe's quantum mechanical approach 

Bethe^ in 1930 developed his theory to sjjfe account 
for the energy loss of fast moving light particles for which 
the classical theory predicted too great an energy-loss. 

Bethe' s'^ treatment is based on the Bom approximation, applied 
to collisions between the heavy particle and the atomic elec- 
trons. Ihe Born approximation requires that 

2 

1 ' ...(3.9) 

where ze and V are the charge and velocity of the moving 
particle. The above condition is well satisfied for hii^ 
velocity and low-charge of the moving particle. The energy 
loss per unit path length is given by 

H ^T dCT/ ...(3,l0) 

where d«ri3 the Rutherford differential cross-section , T is 
. the energy lost and n is the nurrtoer of atoms per unit vol-ume 
of the medium of atomic number Z 2 . Integrating Eq. (S.lO) 
between the limits T^^ and T^^^, 
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dE 

dx 


2 4 
4-rrz e 

Tr2 

mV 


nZ2ln 


T 

max 


T . 
mxn 


. . . (3.11) 


where T is given by the maximum energy transfer by the 
max 3- ^ 

. 2 
incident particle to the electron i.e. 2mV , while T . is 
^ min 

the energy necessary to excite the electrons of the stopping 
atoms. These considerations/ when applied to Sq. ( 3 . 11) / gives 


dS 

dx 


2 4 
4^2 e 

2 

mV^ 


nZ^ 


In 


2 mV^ 

— / 
I 


. . . (3.12) 


where I is the mean excitation potential of the atoms of the 
medium. For relativistic velocities of the incident particle 
the energy loss is given by 


dx 


2 4' 
47rz -e 

m V 





where V/c and c is the velocity of light. 


. . . (3.13) 


III. 2. 4 Bloch' s relation 

Bloch^^ also derived an expression for the electronic 
energy loss rate/ which incorporates the quantum-mechanical 
equation of Bethe, and the classical formula of Bohr, as limi- 
ting cases for small and large values of the parameter zeV('^^V) 
This ' equation has the follov;ing form 


4mz^e‘^ 

mV^ 



+ '^(1) - R “^(l+i 
I 



-RV 


. . . (3. 14) 


dE 

dx 


I 
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where \i/ is the logarithmic derivative of the gamma-f unction 
and RY denotes the ■ real part ofy. For zeVC'RV)^^ 1. Eq. 

(3.14) reduces^ to Eq.(3.7) and for zeV('WV)«^ 1 it is Bethe's 
relation- i.e* Sq.(3.l2) . But for intermediate values of 
ze / (hV) , Eq. (3.14) is approximately equal to either of the 
Eqs. (3.7 and (3.12) . 

Bloch^^ also showed, on the basis of the Thomas-Permi 
statistical approach, that the mean excitation potential I* of an 
atom should be proportional to the atomic number 

l^KZ^ ...(3. 15) 

The enpirically determined values of the "Bloch Constant" K 
do turn out to be approximately the saire for large values -of 
. But for media of low atomic nvuribers, values of K are found 
to be -unexpectedly higher than those predicted by Eq. (3.15) . 

111.2,5 Bohr's semi-quantum mechanical treatment of 
electronic collisions ^ 

Three decades after the development of his classical 
theory of energy-loss, Bohr^ turned to a more detailed consi- 
deration of this problem. Bohr^ dealt with electronic colli- 
sions in detail by taking care of both excitation and ioniza- 
tion as well as -the quantum mechanical "resonance effect" . 

The classical approach in 'terms of a definite irrp act 
parameter is valid only when 1 whereas 7^. '(-C 1 is required 
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for the Bom approximation on which Bethe’ s theory is based, 
where X is given by 

%= 2 zeV(il7), ... (3.16) 

and V is the velocity of the ion of charge ze . The classi- 
cal energy loss Eq. (3.7) v/as put by Bohr in the form 


AE = n AX* Z In (k 


. . { 3 . 17) 


where k = 1.123, is the maximum energy transfer by the 
moving particle to an electron of the medium v/hile is the 
energy transfer in a free collision with inpact parameter 
p = d . The other symbols are defined as follov/s J 

Bg = 2 ■Jrz^e^/CmV^), ...(3.18) 



..'.(3.19) 


* 2 

where tlie ionization potential I = h and V = e /iT- On 

s s o 

the other hand Bethe' s quantum mechanical treatment gives 

T 

A E = 2n AX* £ In ttS , ... (3.20) 


where T = 2m V^ and I is the excitation potential’ of the 
m s 

gth electron. For T yiy I , which is assumed in both the 

m ' s 
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treatments. Sq.(3.19) can be rearranged to give 


'^m 5 2 m 


2z ^ 


, ... (3.21) 


where (2m V^)and(2z ~)have been replaced by and'X.. Eq. (3.21) 
shows that, apart from the factor 3c which m.ay be taken as 
unity in Eq,(3.17), Sqs.(3.17) and (3..-20) become identical 
for %= 1, on the assurrption that the mean ionization poten- 
tial and the mean excitation potential are almost identical. 
For highly charged heavy particles like fission fragments, 
classical treatment, Eq. (3.17), can be utilized to obtain the 
energy loss, which is based on the fact that and b^(^a , 


with 


. . . (3.22) 


and a 


s mU 


. (3.23) 


where m and e are the electronic mass and charge, U is the 

s 

orbital velocity of the s’*'^'^ electron with orbital extension 

a and b is the collision diameter, b <(4 a is always fulfilled 

s ^ s 

for fast moving light particles, but for heavy ions, for some 
values of V, b may exceed a and in which case d , the 

S S , ■ ■ ; , ■ 

inpact parameter p = d^ corresponding to energy transfer Pg' 

will not give an effective adiabatic limit. In such cases 

i h d , where i is the imoact parameter corresponding to 
S / s s . 
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energy transfer I and the energy loss may be somewhat larger 
than given by Sq. (3.17) . Two modes of energy transfer may 
be considered to take into accoimt the collisions for which 
a_^)>))b and a <b. 

(a) Free collisions : If ‘the distance of closest approach 

is of the order of atomic dimensions/ the interaction involves 
the passing particle and one of the atomic electronic and the 
energy transfer lies between and the maximum energy trans- 
fer T^. 'This process is often described as free collisions. 

(b) Resonance effects : If the distance of closest approach 

is large compared to the dimension of the atom, the atom reacts 

as a whole to the variable field set up by the passing particle. 

The result is that the probability of excitation or ionization 

in a single collision is small and the energy transfer T lies 

between D and I (D u^T/I )/ D beina the energy transfer 
s s s^^s s 

corresponding to an iiTpact parameter d . 

s 

Thus the total energy-loss is due to free collisions 
and the contribution from the resonance effects. Two cases 
arise depending on the values of 73 , defined by Eq. (3.16) . 

In the following sections (a) and (b) contributions of free 

collisions and resonance effects v/ill be derived for two 

3 ' " 

values of /pc , as . shown by Bohr. 
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(i) » 1 

For the irrpact parameter p = the corresponding energy 
transfer T = for free collisions, is given by 


I = 
s 


2 ^ 

2z^e* 1 


mV^ i ^ 
s 


. . » (3.24) 


xvhere ze is the charge of the ion of velocity V. Eg. (3.24) 

may be rearranged by using I = 4 mU ^ to obtain 

s 2 s 


- mU ^ ^ 1 ^ 

2 . s ^ Vh' * 2m' ■ 2 


. 2 /2ze^>2 , !r ,2 

^s = ^ 

s 


...(3.25) 


220 *^ 'fl 

Buf%= and ag= / hence Sq.(3.25) reduces to 

‘s' 


ig = 7Ca^ ... (3.26) 

The maximum value of the irrpact parameter p is d correspond- 

.'■''S’- 

ing to energy transfex" in a free collision and is given by 

the time of interaction T of the ion with the electron, v/hich 

should be shorter than the natural period of vibration dd of 

S ' 

the electron. Hie condition for this is 


p = d^ = 


V 

COs 


V h 
I 


or dg= V-n/(|mU^^) 


(2V) 

i ^mU ^ 
s s 


d = ^ a , 
s ‘s s 



. . .(3.27) 
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where = 

• s U 


2V _ , K 


and a. = „ . 

„ s mU 

s s 

d_ 


From Eqs . (3.26) and (3.27) -- is given by 

s 


d 

r = 

s ® 


. . ( 3 . 28) 


From the above equation for, r— (^1, i is always less than d 

% ® s 

and collisions with impact parameter p /i may be considered 

^ , s 

as free and the energy transfer T corresponding to the iirpact 
parameter p is always greater than The energy loss in 

free collisions is given by 

T 

( A E) ^ = n AS^Bg In . . . ( 3 . 29) 

s 

From Eq.(3.26) for’)t^l, ig ) and the collisions with p^ i^ 
are of resonance type. In this case the energy transfer T is 
less than but exceeds corresponding to the inpact para- 

meter d as given by Eq. (3.27) , The contribution of the reso- 

. s 

nance effect to the energy loss becomes 

( AE) . =n AxBg In , . . . (3 . 30> ) 

XT''. J-'' 

■ ; ^ s: ; 

(T /I ) and (I /d ) in Eqs. (3.29) and (3.30) may be obtained 
m s s s 

as follows: 



(2V2 

'u ' 

S 



. . . ( 3.3 i )' 
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1 „ 2 
2 




(— Ar-) ^ ( — ) d ^ 

' ' J-r /vN^ TT ' C 


2z 


Vo' ' -fi/mU' "s 
IT s 


yi~2 _.-l 2 


a- d 
s s 


. . (3.32) 


Using the value of d from Eq. (3.27) in Sq.(3.32)^ I /D 

o s s 

becomes 


s -v. 2 -2 

d“ a. 


. . . (3.33) 


Substituting the values of T /I and I /D in the expressions 

s s s 

for (AE) and ( AE) 


(AS)^ = n Ax In 
( Ae) = n Ax Be In ^ )L~^ 


. . . (3.34) 
. . . (3.35) 


From the above equations it is clear that for 1 ^ ( Ae) ^'^ ( Ae) . 

OC / 

The assumption on V7hich Eq. (3.34) rests is s. 1 and this 

leads to the inequality ■ 


V 


> 


u 




. . .(3.36) 


where 


2V 

U 


. But there are certain velocity regions 


where (^ir) exceed unity and in this case the energy 

loss given by Eq. (3.34) will not be valid. 

from Eq. (3.28.) a i exceeds d_ and the 
values of the inpact parameter p = d for which collisions 
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were "free" no longer extend to i^, since for inpact para- 
meters p>d the duration of the encounter T exceeds , 
s ‘ Co 

3 * S 

Bohr'" showed that the limiting value d of the inpact para- 
meter for which the probability of ionization is still unity, 
is given by 

(d* ) ^ ... (3.37) 

o S 

Using the values of i =>.a and d = >| a , from Sqs, (3.26) 

s s s s s 

and (3.27) , in the above equation. 




. . . (3.38) 


If D * is the energy transf er corresponding to the inpact 

parameter d * , then the energy loss in free collisions is 

s 

given by 


( ZiE) ^ = n Ax in 


m 


.. . (3.39) 


T 

The value of (■^^) may be obtained as follows: 
s 


m 


^s* 2z^e'^/(mV^ 


2 mV ^ .^2 -) (_ i _)2 (^)2 ^*2 

Uz,#J ^ A 

V 


Putting the values of d * f rom Eq. (3.38) and X- 2z „ , 

and 71 =~.in the above equation, becomes 

V. TT IT JJ 


s mU 


m 

D* 


[iV- 


. . . ( 3 .40) 
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and the energy loss is given by 


(AE)^ = nAx Bfc In ^ 1“^ , ...( 3 . 41 ) 

Eq. (3.41) holds for D* and ^ ^ • 

The total energy loss due to free collisions and the resonance 
effects will be the sum of 2qs. (3.34), (3.35) and (3.41), 


(Ae) = n Ax Be l^ln 




where the quantity with in the square bracket is to be replaced 

■% y 

by unity if <^1. 


(ii) %<. 1 

For 1, classical mechanical ideas fail coirpletely 
in accounting for the individual collision effects. But the 
division between the free collisions and resonance effects is 
still valid for impact parameters p^a^ and p^a^/ where a^ is 
the orbital extension of the s‘*^ti electron. The free collision 
part is still given by Eq. (3.29) \-jhlle for p/a , which corres- 
ponds to distant collisions (resonance ef fee ts) , the energy 
loss is 

A ' 

( AE) = n Ax Be In. , ... (3.43) 

^ s 

where A ' is the energy transfer corresponding to the inpact 
s 



32 


parameter p = a and is given by 

5 


A “ 
s 

2 z^e^ 

1 

mv 2 

a 2 
s 

7v * 

or A 

s 

II 

• 

H 

s ' 

A’ 

^ s 

and — , 

s 

I 

_ s 

D ' 
s 

%2 


2 zV 2 A-t o 1 

( O) (JL) 2 _1_ 

V '•mU ^ 2 

s a 


2 “V 


(3.44) 


. . . (3.45) 


where is given by Eq. (3.33) . Incorporation of Eq.(3.45) 


in Eq.(3.43) leads to 


( E) = n A X In .2 ... ( 3 . 46 ) 

s 

. and (AE)^=nAx Bcln “^2 (3.47) 

X s 

Thus for%(^l/ the total energy loss is given by 

( AB) = 2 n Ax B 4 In ...(3.48) 

• s 

Eqs. (3.42) and (3.48) combined together give a conprehensive 
energy loss equation 

Bfe (S In ^ g ^ ' .. . (3.49) 

s 

where the terms within the square brackets, if less than tinity, 
should be replaced by unity and the . summation of the logarith- 
mic terms include the velocities of all those orbital electrons 
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which interact with the roving ion at a given ion velocity V. 
The other symbols are defined as follows: 


Be 

= 2Jfz^ej^(mV^) 

. . . ( 3 . 5o) 

s 

= 2V/U 

S 

...(3.51) 

X- 

= 22V / V 
o 

. . . (3.52) 


where m and e are the electronic mass and charge respectively/ 
^s orbital velocity of the s^^ electron in the atom of 

the medium/ = e^/lT , ze is the charge of the ion and "h, is 
the nurriber of atoms of the medium per mit volume. 


111.2.6 LSS Treatment 

For the energy loss of charged particles in stopping 
media with atomic numbers less than 10/ Lindhard et al.^^ 
worked out a theory based on the Thomas-Ferrni statistical 
model pf the atom. They expressed the energy E of the moving 
ion and the range R in terms of dimensionless parameters £ 
and P defined by 


€= 




B a 


^’2) 


and P = R n 


4 ^a^ Ml 


TM 4- M 


s 2 


( 3.53) 

(3.54) 


vjhere and , Z^^ are the mass nioiribers 


and the atomic 


iri 


34 


numbers of the moving ion and the stopping medium respectively, 
n is the number of atoms of the medium per unit volume, and 
a is the screening parameter and its value is given by 

a = a^ 1^0.8853 

Sq being the Bohr radius of the hydrogen atom. The electronic 
stopping is assumed to be proportional to €. 



.. . ( 3 . 55 ) 
. . . ( 3 . 56 ) 


III. 3 NUCLEAR COLLISIONS . 

The energy loss of heavy ions in elastic collisions 
with the atoms of the stopping medium tabes place only near 
the end of the ion range, where the probability of electronic 
excitation becomes small. Bohr^ showed that such two-body 
'nuclear' collisions may be described classically if the 
collision diameter of the moving ion greatly exceeds the 
de Broglie wavelength. This condition is given by 




57) 
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where and Z^ are the atomic numbers of the moving and 
struck particles respectively and V is the ion velocity. 


To describe a 'nuclear' collision between two heavy 

. 18.W 

ions, the Rutherford scattering^can be applied directly. In 
terms of energy transfer the differential cross section 
for collision between two particles is given by 


d(r(S/ T) 


2 2 4 

7rz/z_^e^ M. 

± ^ (-i) 

E 

2 T 


2 ' 


« • * { 3 • 58 ^ 


where and are the masses of the moving and struck 
particles. In the derivation of Eq. (3.58) it is assumed 
that the collision is occurring between two point charges with 
central fields of force. In the case of collisions between a 
heavy ion or a fission fragment and an atom of the stopping 
medium the presence of bound electrons on each nucleus shields 
the nucleus and in that case the simple Coulomb potential has 
to be modified to take care of the screening effect of the 
bound electrons. One such potential V(r) is given by 

V(r) = (Z^Z^eVr) exp(- r/a) 

where a is a screening parameter and r is the distance. The 
applicability of the Rutherford law depends critically on 
the ratio between ‘ a', and the collision diamoter ' b' for an 
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unscreened field -^ = or b and also 7 C)'^I, 

then the enounter would be a very intimate one and the inter- 
action would be confined mostly in the unscreened part of the 
field. For strong, screening .^^^1, a is small so that colli- 
sions may be considered as occurring between two rigid spheres. 
If R(S) is the hard sphere radii which depend on the energy,, 
the differential scattering cross section is then given by 


d<r{E, T) = 4 '!r[R(E)J ^ . . . ( 3 .60 ) 


where T 


max 


4 


(M^ + M^) 


E 


The differential cross section is related to the rate of energy 
loss by 


dx 


n 


1 


max 


T dcr<E/ T) . . . (3 .61) 


min 


By a simple geometric treatment of the scattering process/ 
Bohr^ showed that 


mxn 


e"^ M n 
1 2 /__iN jL 

2 M/ E f 

o Z 


. . . { 3 .62 ) 


I' 
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and the rate of energy loss due to nuclear collisions is 
given by 


dx 


2 2 4 

2iT2/z/e 

= ( r ) n In 


^ 1^^2 


^r2 

V a 


ZlZ2e^ 


( 3 . 63 ) 


where a is the screening parameter and all the other symbols 

o 

have been defined earlier. 



CHAPTER IV 


DSVELOI^IENT OF RANGE-VELOCITY 
AND STOPPING- POWER EQUATIONS 
FOR FISSION-FRAGMENTS AND PARTI- 
ALLY STRIPPED HEAVY IONS 


IV. 1 DERIVATION OF THE RANGE-VELOCITY RELATION 


The energy- loss per unit path length of a moving ion 
of charge ze and velocity V is given by Eq. (3.49) . 


dx 


O TT ^ ^ 

2^z e n 


mV' 


inin 


U 


wi 


dn(U ) + 
s 


u 


Ktsr}"' 


... (4.1) 

The above equation differs from Eq. (3 .49) only in the replace- 
ment of the summations by integrations , The integration should 
include the orbital velocities of only those electrons which 
interact with the moving ion at any given ion velocity V. 
Because of the continuous capture and loss of electrons by the 

decelerating ion its charge fluctuates and in that case ze in 

■ '0f f 

Eq. (4.1) has to be replaced by an 'effective charge' Z of 
the ion. The evaluation of Z and the integration of the 
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logarithmic terms were done by Bohr'' in the case of a "heavy" 
ion passing through a "heavy" medium. 


IV . 1 * 1 Heavy solid media 

Two assumptions were made by Bohr;^ 

(a) The number of orbital electrons n(u ), whose velo- 

s 

cities are less than a given velocity U , is given by 

s 


n(u^) 


u 



. . (4.2) 


v/here v is the "effective" quantum number for the outer 

o 

orbital electrons. In the case of heavy atoms (z'^-40) the most 
firmly bound electrons, belonging to K, L, etc shells, move in 
a field which is approximately coulombic and have values of 

V very nearly equal to 1,2,... etc, respectively. Further- 

s 

more the most loosly bound electrons for which n(U )~1 again 

S' 

have values of of the order of unity. But over a large 

intermediate region, v has a flat maximum corresponding to 

s 

1/3 . 

values close to Z , where Z is the atomic number of the atom. 

This result is in conformity with the analysis of electron 

binding by the Thomas-Fermi statistical model, n(U ) was 

3 

expressed by Bohr as 

n(U ) = , . ..(4,3) 

o 

holding for For values of outside this region 

Eq, (4 . 3) is no lon^^ valid and in that case n(U^) can be 
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represented by some function of Z and U ^vrfiich for U v is 

^ s s o 

of the order of vinity, while it approaches Z for U — ZV . 

s o 

(b) The second assumption^ made by Bohr, is that the 

number of electrons lost by a heavy ion at a velocity V is 

equal to the number of orbital electrons with velocity less 

0f f 

than V . Ihus, Z can be given directly from Sq. (4 . 3) . 


= n(U = V) = Z^^^ — . 
s V 

o 




The number Z ' was estimated for the heavy fission fragments 
. . 30 

slowing down in gases. The effective charge of the frag^ 
ments in solids is difficult to evaluate, but it is expected to 


be higher because of much shorter time between collisions, 

31 

compared to the time of de- excitation. Bohr and Lindhard 

0ff . 

made a qualitative estimate of Z which is given by 



where k = 1.5 in solid and k = 1 in gaseoios media. The evalu- 
ation of two integrals of Eq. (4.1), was done with the help 
of Eq. (4.3) . Ihe lower limit was taken as U^= O and the 
upper limit U ‘ corresponds to those values of U which make 
the logarithmic ter^Bzero. Considering these two integrals 
separately, one obtains 
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J, = 


U =U' 

' s s 

3 in^ 


U =0 
s 




dn(u ) 

s 


z V3 m 

fa 

V 


j ln[h2[xfjda 


= 4 ~ yT^, 


and 


U =U ' 
s s 


J 


In 


Us=0 


{vitn 


V 


dn(U ) 
s 


. . . (4.6) 


. (4.7) 


To evaluate the above integral the integration has to be 

performed in two parts from 0 to 2 V%^ for ^ <^1 and 2 VtC^ 

^s ifl 3 Using these considerations in Bq. (4.7)^ 

s 2 

becomes 


J„ = 


r .2VX"^ 


2 V 


L 


InV^t ^fdU + 
s \ s 


^-1^" [t ] ' '5"s] 

r -b ^ S 


r- 


V 


1^4 4 In 


2vyL 

vl/3 


4 (ex 


- 6%^ - 4% 




— i — - [3>^^ - tC^] 

r\ -* 


. . . (4.8) 


The stum of the integrals in Eq. (4.1) is given by 

tc^] 


J = = 


2Z 


V 


. (4.9) 


Use of Eq. (4.9) in Eq. (4.1) leads to the following energy loss 
equation: > 
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dx 


4Tre^N 

-5 

mV “ 


o 





,. (4.I0) 


Using Eg. (4»lo) and dE — A^m^VdV, whene is the fragment 
mass number and m^ is the nucleonic mass and substituting the 
values of all the constants, one obtains the range- velocity- 
equation 

A A dV 

p(3x = ^ ■■ 

127.3 X 10^^ Z2^'^^j4.7622(kZ^^/^)5/3 + 

...(4.11) 

where A^ and Z^ are the mass number and atomic number of the 
moving ion while P , A^ and Z^ stand respectively, for the 
density, mass nxuriber and atomic number of "the mediiam. The 
total range or "integrated range" might be obtained if one 
could integrate the right-hand side of Eq. (4,11) down to zero 
velocity. However, Eq.(4,ll) is not strictly valid at , 

since at these velocities the ion is practically neutral and 

stopping will be mainly due to screened field type of interac- 

12 ■ 
tions. Niday advanced the plausible argxjment that as far as the 

projected range' in the direction of the moving ion is concerned, 

which corresponds to the experimentally measured range, it is 

sufficient to evaluate the integral between the initial ion 

velocity and the velocity at which -fche ion becomes almost 

neutral. At there would be insignificant addition to the 

projected range because of the very large angle scatterings 
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associated with nuclear collisions. Hence, the experimentally 
determined range may be related to the initial velocity by 
the following expression: 


0 

R(mc[/cm ) = j Pdx = 


^1^2 V 


127.3 Z2^^^j4-7622(kZ^^/^) 


. - . (4.12) / 


where both and are to be expressed in units of lo^ ciry'sec. 
and all other symbols have been defined earlier, Niday^^ showed 
that with k « 1, the ranges of various fission products in 
uranium can be calculated quite accurately. Since k 1 is 
different from what Bohr and Lindhard suggested, all availa- 
ble range data for fission products were analyzed in heavy 
solid media in order to investigate if there is any systematic, 
medium-dependent variation in the value of k. The procedure 
used for calculating the value of k is as follows. For a given 
fission product of mass number of knov/n range R, the initial 
velocity is first calculated. For this purpose, the pre- 
cursor fragment mass nxamber A^' is first determined by means 
of the experimentally measured prompt-neutron number as a 
fxanction of fragment mass number. From the experimentally 
measured fission fragment kinetic energies, existing in 
literature and A^, ' , the values of have been calculated. 

Since the prompt neutrons are emitted almost isotropically 
fiXDm the moving fragments, the velocity of the fission product 
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of mass nuxriber is identical with the velocity of the pre- 
cursor fragment of mass nioiriber A^'-. Ihe atomic number 
(also termed the nuclear charge) of the fission product of 
mass number A^ is given by the "most-probable charge" Z»(A.) 
which has been calculated from Mukherj i' s^^ prescription 

,A -A ’ 

for light products = Z^- ...(4.13) 

A ^ 

for heavy products -.-(4.14) 

where A and Z are the mass and the atomic numbers of - the 

* # 

parent compound fissile nucleus. Using appropriate values of 

Vf/ and, A^ in Eq,(4.l2) the values of k for different 

fission products in a particular mediuni have been obtained. 

12 

The range values have been taken from Niday for uranium, 
Alexander and Gazdik^^ for Au, Almodovar et al.^^ and 

I 34 35 

Hontzeas and Blo(|ik for tungsten and Smith and Prank for 
zirconium. In all cases, the precursor fragment mass number 
A^‘ has been obtained with the prompt neutron emission data 

O 

of Apalin et al. The fragment kinetic energies have been 

37 

taken from Schmitt et al. Table IV. 1 shows the k valiJes 
obtained with different fission products in uranium. For 
other heavy catchers like Zr, W and Au/k values are listed 
in Table IV. 2. The arithmetic mean values of k obtained from 
Tables IV i and IV. 2 in U, Zf, W and Au together with the 
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TABLE IV. I 

Values of k calculated from Eq. (4.12) for 
different fission-products in uranium 


Nuclides 

Velocity 
(in linits of 

8 

lo cm/ sec) 

Most probable 
charge (Zp) 

k 

pmnupmigiiiiiiiiiiiiiiip 





14.9407 

37.00 

1.0950 

®"sr 

14.6398 

35.46 

1,0544 

^°sr 

14.5750 

35.95 

1.0359 

Sr-Y 

14.4956 

36.33 

1.0551 

93^ 

14 . 3405 

37.11 

1.0644 

95„ 

Zr 

14.1478 

37.92 

1 .0602 


13.9922 

38.75 

1.0585 

^%Io 

13.8829 

39.54 

1.0700 


13.5270 

41.20 

1.0562 


12.9561 

42.42 

1.0541 


12.3880 

43.68 

1.0611 

111. 

Ag 

11.9988 

44.45 

1.0819 

112 

Pd 

11.8298 

44.98 

. 1.0826 

115^, 

Cd 

11.3200 

46.09 

1.0610 

■125 

^^Sn 

10.9786 

48.49 

1.0897 

127 

'^'Sb-Te 

10.9596 

49.09 

1.0607 

129^3^ 

10.9296 

49.98 

1.0524 

132^^ 

10.6825 

51.21 

1.0492 


9.9582 

53.34 

1.0347 




contd. 


TaJple IV. 1 (contd.) 
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1 

2 


9.6258 


9.4863 


9.1912 

144 

Ce 

9.0879 

147-^ 

Nd 

8.9047 


4 


54.69 

1.0457 

55 .08 

1.0498 

56.04 

1.0383 

56.43 

1.0346 

56,43 

1 .0307 


TABLE IV. 2 


Values of k calculated from Eq.(4.12) for 
different fission products in Zr, W, and Au 


Nuclides 

‘Wzr 

‘«W 

S^Sr 

0.9739 

0.9587 

^"zr 

- 

0.9637 

■ 99.. 

Mo 

1.0040 

0.9962 

“^Ag 

- 

0.9933 



0.9858 

131 j 


0.9818 

I32^e 

“ 

i.oolo - 

140„ 

Ba 

0.9876 

1.0030 0.9637 


TABLE 

IV. 3 

Mean values 

of k in different heavy media 

Medium 

Atomic number k 

‘V 

Mean 

Zirconium 

40 

0.98840.05 

Tungsten 

74 

1.0004<).003 

Gold 

79 

0.98010.03 

Uranixm 

92 

1^06 40.03 
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maximum deviation of. an individual value from the mean value 
are .shown in Table IV. 3. Since there are finite errors asso- 
ciated with the measured ranges, it is possibly justifiable to 
assume that in heavy solid media and possibly in all solid 
media k may be considered as unity. This analysis supports 

Bohr' s original assumption that ^ , independent 

o 

of the medium being penetrated. Hence the ranges in all heavy 
media can be obtained from the following expression: 


R(mg/ cm^) 




127.3 Z. 


^^.7622 




. . . (4.15) 


IV . 1 . 2 Light solid media 

In the case of light solid media, the main problem lies 
in the integration of the logarithmic terms in Eq. (4 . 1) .Since 
the minimum value of in the case of fission products is 
close to 40, the use of the Thomas -Fermi approach is justifi- 
able and Eq.(4.4) can be used to calculate the effective 
charges. Hov/ever, for media like Be, C, A1 and Si , 

Eq.(4.3) is certainly not api^lic able, since no statistical 
approach can be valid for systems containing only a few elec- 
trons. Bohr^ used Eq. (4. 3) to evaluate the summation terms 
of Eq. (3.49) for heavy media. In our analysis, we have procee- 
ded with the assumption that in the case of light atoms (Z^40) 

take n(U ), the number of electrons having yelocity 
s 


one may 
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less than a given velocity U to be 

s 

U 

n(U^) = f(Z) ^ , ...(4.16) 

o 

where f (Z) is a function of Z whose forr has to be determined. 
With this asscim|3tion and following Bohr* s^ procedure of repla~ 
cing the summation terms in Eg. (3.49) with integrals, and 
making appropriate provision for the non-participation of many 
of the orbital electrons of the light media, one gets an expre- 
ssion analogous to Eg. (4.9) for the integral J. 


J = 2f(Z ) + >1] ...(4.17) 

o ^ 


After siibstitution of the numerical values of the physical 
constants and the replacement of the integrals in Eg, (4.1) by 
Eg. (4.17), the energy loss eguation becomes 


M = 1.327 [4.7622 Z + z/Zq V 

dx mg ^2 ^ ^ ^ ' 


.... ( 4 . 18 ) - 


12 

Further using dE =AVdV and Niday' s assumption for the 
lower cut-off velocity V^, the range-velocity eguation becomes 


R (mg/cm ) = 


^1^2 ^^i-V 


127.3 [4.7622 Zf 


5/9 


+ Z ^ " 

^ r 


. . (4.19) 
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In order to evaluate fiz^), the san^ procedure has been 
followed as in the case of heavy media for evaluating k. 

Range data for various fission products in' a particular 
medium along with the corresponding values of Vj_ and 
were put in Eg, (4.19) to obtain individual values of fCz^) in 
one medium. Ranges in A1 were taken from Alexander and 

Gazdi'k, Aras et al.. Brown and Oliver, and Nakahara 

39 1 Q 

et al., in carbon from Chinaglia et al,, and in Be and 

• 4-0 

Si from Demichelis et al. Only one datum on range in Cu 

17 

is available from Segre' and Wiegand for the average fission 

O -DC 

product from the thermal neutron fission of U. since the 
fission products in equal and maximum yields have mass and 
atomic numbers 140, 99 and 54.70, 39.44 respectively, the 
values of A^ and Z^ for the average fission product may be 
taken as 119.5 and 47.07 respectively. Using similarly an 
average velocity '7^, £(^ 2 ) for Cu has been calculated from 
Eq. (4.19) . Table IV. 4 shows the values of ^(S^) obtained for 
different fission products in Be, C, Al and Si, The arithmetic 
mean value of £(^ 2 ) together wi-th its maximum deviation from 
an individual value for any particular medixm is shown in 
Table IV. 5. Zirconium has also been included among the 
light media to locate the dividing line betv/een the ••light*’ 
and 'the "hea'vy" media. The plot of f (Z 2 ) "vs Z^ , shown in 

Fig. 1, shows a slope of 0,28 and thus f(Z 2 ) may be very 
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TABLE IV. 4 


£( 22 ) Values for Be, C, Al, and Si obtained from 
Eq. (4.19) using different fission products 


Nuclides 

Be' 

C 

A1 

Si 

8’Rb 

- 

- 

1.6019 

- 

S^Sr 


- 

1.5390 



- 

- 

1.5190 

- 

92y 

- 


1.5775 

- 

952^ 

0.7412 

0.9110 

- 

1.6153 

97 

Zr 


- 

1.5728 

- 

^^Mo 

0.7149 

0.9007 

1.5815 

1.6035 


0.7185 

0*9546 

- 

1.5794 

'“Ag 

\ 

- ■ 

1.5937 

- 


- 

- 

1.5850 

- 

isij 

0.7016 

0.9103 

1.5709 

1*6008 

132tq_j 

0.6823 

0.8660 

- 

1.5768 


- 

- 

1.6190 

- 

l^^Ba 

0.7221 

0.8502 

1.6088 

1.6148 

I4lce 

0.6953 

- 

- 

- 


- 

,0.8253 

- 

1.5911 


0.7043. 

0.8128 

- . 

' -■ 
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table IV. 5 

Mean f (Z^) values for different light media 


Medium 


Atomic ntunber 


f(Zj) 


Hydrogen 

1 

0.275 + 0.012 

Helium 

2 

0.347 + 0.009 

Beryllium 

4 

0.710 + 0.03 

Carbon 

6 

0.879 + 0.07 

Nitrogen 

7 

1.080 + 0.02 

Neon 

10 

1,166 + 0.02 

Aliominium 

13 

1.582 + 0.02 

Silicon 

14 

1.597 + 0.02 

Argon 

18 

2.182 + 0.04 

Copper 

29 

2.64 

Zirconium 

40 

3.23 + 0.33 


1.1. T. K^NFUR 
QENTRA^ library 

Acc. N«. i ■' 


t (Z 2 ) 



0 1 2 3 4 5 6 7 8 9 10 11 12 


FIG. 1. Plot of/(2^2) VS Zy^ for various media. 
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accurately given by 


f(Z2) = -0.28 


. . (4.20) 


Further, from n(U^) = ^ = O.SSZ ^/^ ^ , it is found 

_ o o 

that = 45.5 constitutes the dividing- line, between the light 

and the heavy media. Ihus/ the number of electrons n(U ) 

s 

having velocity less than a given velocity U in the case of 

s 

media with Z<^45.5 is given by 


n(U^) = 0.28 


2/3 ^s 


...(4.21) 


Similarly, the effective charge of an ion of velocity 
V(Z^5.5) is given by 


~eff 2/ 3 V 

Z^ ^ = 0.28 Z^/ — 


. . . (4.22) 


Eqs . 4.15, 4.19 and 4.21 lead to the following comprehensive 
range-velocity equation in solid media: 


R(mg/cm'^) = 




127.3 f(Z2) [4,7622 ^f(Z^)j^/^ + f(Z^)] 


, . . (4.23) 


where A^, and Z^, Z^ are the mass numbers and the atomic 
numbers of the fission product (or heavy ion) and the medium 
respectively, is the initial velocity of the moving ion 

and V = 2 .18x10® cm/sec . Both V . and V are to be expressed 

o ' . ■ " ^ 

innnits of 10® crr/sec.> f(Z^) or f(Z 2 ) or in general f(Z) is 
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given by 

f(Z) = 0.28 Z 445.5 ...(4.24) 

f(Z) = for Z^45.5 ...( 4 . 2 s) 

IV. 1.3 Light gaseous' media 

Range data for fission products from the fission of 

2 33 . 

U in the light gases H 2 /He, Ne and Ar have been taken 

4 1 

from Pertzhak et al. and the values were calculated 

using Eg. (4.19) . The initial velocity and the precursor 
fragment mass nxjirber have been calculated with the data of 
Pleasonton"^^ and Apalin et al.^^. Table IV. 6 shows the results 
of the calculation for £( 22 ) for different gases. 

TABLE IV. 6 


Values of f (^ 2 ) for gases obtained from Eg, (4.19) 


Nuclides 

Hydrogen 

Helium 

Nitrogen 

Neon 

Argon 

91^ 

Sr 

0.2788 

0.3545 

1.0791 

1.1635 

2 . 1478 

CM 

o\ 

0.2784 

0.3558 

1.1075 

1.1556 

2.1619 

S’zr 

0.2864 

- 0.3512 

1.0986 

1.1790 

2.2060 


0.2740 

0 . 3485 

- 

- 

2.2220 


0.2555 

0.3261 

1.0552 

- 

2.1764 
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The arithmetic mean values of f(Z^) thus obtained are sho;vn 
in Table IV. 5. From Fig. 1, it is evident that except for 
Ar/ f (Z^) for other gases ^2' are accurately 

given by Eq, (4.20) . In the case of hydrogen and helium this 
agreement may be fortuitous^ but this would be useful in 
calculating the ranges in compound media like rrylar and 
collodion, which contain a number of hydrogen atoms, as shown 
in Sec. IV. 1.5. 


IV, 1.4 Ranges of accelerated heavy ions 

Although there should not be any difference between a 
fission product and an accelerated heavy ion as far as the 

applicability of Sq.(4,23) is concerned, it would be inheres- - 

■ 22 
ting to verify this. Bridwell and Moak have measured the 

127 79 

ranges of I and Br ions accelerated to different ener- 
gies in the form of R{e)-R(1o), where R(e) and R(10) are the 
ranges corresponding to energies E(MeV) and lo MeV respectively. 
Figure 2 shows these values along with our calculated values 
for the four- elements beryllium, aluminium, nickel and gold. 
Hower and Fairhall"^^ have measured the energy degradation of 
^Be ion during its passage through different thickness of 

aluminium and gold. In the case of gold the calculated value 

2 

of the range R(6.8 MeV)-R(2.2 MeV) is 4.504 mg/ cm , whereas 
the experimental value is found to be 4.50 me/ cm , R(6.8 MeV) 
and R(2.2 MeV) being the ranges corresponding -to initial 



[r(E)-R (10)] ( mg/ cm^ ) 
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FIG. 2. Plot of 
i?(10) against {2E/M)^^’^ for 
'^®Br and ions in Be, Al, 
Ni, and Au. R{E) and il(lO) 
are the ranges correspond- 
ing to initial energies E 
MeV and 10 MeV, respec- 
tively. J and 5, experi- 
mental values from Brid- 
well and Moak, Ref.2i2. ; 

— , calculated values. 


57 


energies of 6,8 MeV and 2.2 MeV. In the case of aluminium 
no such reliahle inrormations could be extracted from the . 
experiment. 

IV-.1.5 Compomd media 


The stopping power equation given by Sq.(4.18) may be 
written as ' , 


^ = 
dx 


a n 


fCz^) 


[4.762 ^f(Z^)j 


5/3 



(4.26) 


where oc is a constant and n is the number of atoms of the 
medi^Im per unit volume and £(2^,) and fCZ^) or f(Z) is given 
by Eqs , (4.24) and (4.25) . If a. compomd is made up of 
number of atoms of the i atomic species per molecule, then 
the actual stopping power would be the sum of the stopping 
powers due to all the individual atoms in each molecule. Hence, 


„ PN 

dx A 

c 


\l .7622[f (Z^)j 


5/3 


+ f (Z^) 


Jv, ...(4. 


27) 


where P is the density of the medium, N is the Avogadro number, 

A is the mass number of the molecule A_ = £ X.A. ; Z. and 
c ^ i 

til. 

A^ being the atomic number and the mass number of the i 
atomic species. Hence, for a compound medium, the stopping- 
pov;er and range-velocity equations are 


dx 


f(^ 

A 


[4.7622[f (Zj^)| + £( 22 )^ V, . . .(4.28) 


1.327 
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R = 


A A (V.-V ) 

1 C 1 


127.3 f(Z^) 1^4.7622 3 + f (Z^) ] 


, . . . (4.29) 


where f(Z ) S f(Z^) and A = E X.A. and f(Z) is defined 
by Eqs. (4.24) and (4.25). The calculation of ranges in 
raylar (C^^HgO^) , collodion (^^2^17^16^3^ and UF^ would 
involve some assumptions regarding f (Z^) . For carbon and 
uranium as also for hydrogen, f(Z^) is known. Since nitrogen 
also conforms to f(Z) = 0.28 , it is assumed that this 

would also be true in the cases of oxygen and fluorine. Using 
the' appropriate mass and atomic numbers and the initial velo- 
city for the average fission product from the thermal neutron 
235 

fission of U as indicated in Sec. IV. 1.1, range in collo- 
dion has been calculated by means of Eg. (4.28) and is shown 

in Table IV. 7, along with the corresponding experimental 

17 

value of Segre' and Wiegand. The calculated and the experi- 
mental ranges for the "average" fission product in aluminixam 
are also shown in the same table as an additional verification 
of the correctness of our procedure for calculating A^, Z^, 
and V^ of the "average" fission product. Table IV. 7 also 
lists the calculated values of the ranges of fission products 

of mass number 99 and 140 in Csl and mylar along with the 

44 

corresponding experimental values of Suzor and Gumming and 
Crespo.® Figure 3 shows the plot of the calculated values of 



Calculated and experimental ranges in mg/cm^ in compoxjnd medi 
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Gurnming 

Suzorl^4 
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0 0,4 0.8 1.2 1.6 2.0 


✓2E/M (MeV/amu)'^^ 

FIG. 3. Plot of i2(£')-i?a0) vs ( 2 J 5 /M )^^2 
ions in UF 4 . R(E) and R(10) are the raises corre- 
sponding to initial energies E MeV and 10 MeV, respec- 
tively. The sjonbols with the error bars show the ex- 
perimental values, while the straight lines represent the 
calculated values ( — •). 


A 
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R(E)~R(l0) in UF^ vs ( 2E/xM) along with the corresponding 
experimental values 2 2 for "^^Br and ions, where R(E) and 

R(10) correspond to ranges at initial ion energy of E and 
lo MeV respectively, ana. M is the mass number of the heavy ion 

IV. 2 STOPPING-POWERS 
12 

Niday' s assumption, which has been used in the deri- 
vation of the range-velocity relation that at V<y^ the pene- 
trating ion does not travel much in its original direction, 
needs an independent verification. The best procedure would 
be to determine whether the stopping powers given by Eq. (4 .1) , 
after the incorporation of the phenomenological deductions 
represented by Egs. (4.16), (4.24) and (4,25), does predict 
the stopping powers in different media at different energies 
which are in agreement with the available experimental data. 
The comprehensive stopping-power equation is 

( !^-ey - -C il£) = 1.327 ~^f4. 762 ^f(Z)f/^ + f(Z)lv .. (4,30) 

dx mg 1_ L i J j 

where V is the velocity of the ion in units of lo cn]/sec at 
which the stopping power is being calculated and ail other 
symbols have been defined earlier. In the case of a compound 

of known molecular formula, the observed stopping power would 
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be the sum of the stopping powers of the individual atoms in 
the molecule and is given by Eg. (4.28). Comparison of the 
experimental and calculated values of the stopping powers are 
shown in the following sections. 

IV. 2.1 Fission fracfments ' 

Figure 4 shows the stopping pov/ers of diffei^ent fission 
products from the spontaneous fission of ^^^Cf calculated with 
Eg. (4.30) along with the corresponding experimental values 
of Bridwell and VJalters^*^ (using thick foils of silver and 

gold) . All calculated values are at V = 1.38 x lo^ crr/sec 

9 

in the case of carbon, V = 1.35 x 10 cn/sec in the case of 

9 

silver and V = 1.32 x io cn/sec in the case of gold. Ihe 

experimental errors, which are of the order of 25-30/ in the 

cases of silver and gold and is somev/hat higher in the case 

9 

of carbon, are not shown in the figure. Kahn and Forgue 
have measured the energy degradation of the median light and 

the median heavy fission products from the spontaneous fission 

/p c* Q 235 

of Cf and the thermal -neutron fission of U, using foils 

of different metals of varying thicknesses. Using these 

authors' values for the mass and atomic numbers of the fission 

products, the experimental stopping powers have been obtained 

in the following manner. From the incident and emergent 

energies E^ and of the fission products before and 



64 


passing through a foil of thickness * t' , the experinjental 
stopping power is obtained as (E.-S^)/t. Ihe corresponding 

calculated from Sq. (4.30) using a value of V corres- 
ponding to the mean energy j (E.+ S^) . Figures 5 and 6 
show these values in aliominium, nickel, silver and gold. In 
a similar manner the energy- loss data of Muller and 
GonnGHwGin ]ia.V 0 Idgsh tirGS-tGcl to obtain tha ©xpGjtriitisntal 
and calculated stopping powers for the median light and 

median neavy fission products from the spontaneous fission 
2^2 

of Cf . Figure 5 includes these results only in the case 
of copper; for the other metallic media the agreement between 
the experimental and calculated values are excellent but are 
not shown in the figure. Only one experimental datum on 
(— ) for the median light fission product from Cf could 
be obtained from the data of Kahn and Forgue^ in the case 
of U-Pd alloy, which is shown in Table IV. 8 . The alloy 
containing 20-v;t,% U and 80-wt.% Pd has the molecular formtila 


U-Pd. 


■ 8.88 


and using Sq. (4 .28) , the stopping power has been 


calculated and is shown in Table IV. 8 . 


IV. 2 . 2 Accelerated heavy ions 

Table IV . 9 shows the experimental values of the stopping- 
powers of accelerated IN ion . as measured by Nakata along 
with the corresponding calculated values with Eg. (4 • 30) in 
alirniniuni, nickel/ silver, and goldv 
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0 4 8 0 4 8 12 


v/FCMeV)'^* 

FIG. 5. Plot of the stopping powers (dEldX) vs /E for 
the median heavy and median light fission products from 
the spontaneous fission of ^®^Cf in Al, Ni, Cu, Ag, and 
Au. The symbols with the error bars show the experi- 
mental values, while the straight lines represent the cal- 
culated values. The experimental values in Al, Ni, Ag, 
and Au are from Kahn and Forgue (Ref. 3 ) and in Cu from 
Muller and Gbnnenwein (Ref. 




/E(MeV}''^ /E(MeV)'''* 


FIG. 6. Plot of the stopping powers (dF/daOvs for 
the median heavy and median light fission products from 
ttie thermal neutron fission of in Al, Ni, Ag, and Au. 
The symbols with the error bars show the experimental 
values taken from Kahn and Forgue (Ref. 5 ) and the 
straight lines represent the calculated values. 
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table IV. Q 

Calculated and experimental values^ ^ of the 
stopping powers of, ion in Al,Ni,Ag and Au 


Stopping 

medium 

Mean 

energy 

(MeV) 

(dE/d ) Cal. 
(MeV cmVmg) 

( dE/ dx) Exp. 
(MeV cm^/mg) 

% Uncertainty 
in exp . ( dE/dx) 

A1 

( 8.00 
( 

5.23 

5.00 

5.59 


( 6.90 

4.85 

4.40 

6.04 

Ni 

( 6.65 
( 

3.59 

. 3.60 

6.28 


( 4.06 

2.80 

2.75 

8.11 

Ag 

5 .50 

2.53 

2.65 

6.85 


( 5.85 

1.69 

1.85 

6.53 

Au 

( 

( 6.75 

1.82 

2.00 

6.16 


Only two values of the stopping-powers in the cases of alumi- 

niiom, nickel and gold are shown in Table IV. 9 , while in the 

case of silver only one datum could be obtained from the, experi' 

dE 

ment. The overall uncertainty in (— ) has been obtained from 
the errors involved in the measurement of the foil thickness 
and the energies of the incident and emergent ions. Shane 
and Seaman^® measured the stopping powers of 18.4 MeV and t 
16.9 MeV ^^Ne ion in aluminium using a new technique in which 
the Doppler shift of Yvrays from a moving ion is used to 
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measure the velocity of the ion. Table iv. lo shows the 
calculated values of the stopping powers with the corres" 
ponding experimental values . 


TABLE ly .ln 

Calculated and experimental values"^^ of the 

20 

stopping powers of Ne ion in aluminium 


Energy 

(dE/dx) Cal. 

( dE/ dx) Exp. 

(Me^ 

(MeV cm2/mg) 

(MeV cmVnig) 

18.4 

9.46 

9.4+0. 6 

16.9 

9.05 

10,6+1.8 


Figures 7 and 8 show the experimental values of the stopping 
powers of accelerated I and Br ions as measured by Mo ah 
and Brown*^"^ at various energies along ■with ■the corresponding 
calculated values with Sg. (4,30) in beryllium^ carnon/ al'umi™ 
nium, nickel, silver and gold. In Fig. 9 the values of the' : 

stopping pov/ers of andl^Br ions in UF^ are shown along 

■ ■ ' ■■ ' ' 22 ' 

with th§ corresponding measured values of Bridwell and Moak. 

Figure 10 represents the experimental values of the stopping 
powers of aoceleratefl ”®u Ion in carbon, aluminlnm, nickel, 
silver and gold from Brown and Moak^® and the corresponding 
values calculated with Eg. C4.30) . Ihe experimental values 
have +10% error which are not shown m the figure. 
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FIG. 7. Plot of the stop- 
ping powers {dEfdX) of 
ions vs 'Te in Be, C, Al, 

Ni, Ag, and Au. The ex- 
perimental values from 
Moak at\d Brown (Ref. 47 ^ 
are shown by the symbols 
with the error bars, while 
the straight lines represent 
the calculated values. 
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FIG. 8. Plot of the stop- 
ping powers (dE/dX) of 
ions vs /E in Be, C, Al, 

Ni, Ag, and Au. The ex- 
perimental values from 
Moak and Brown (Ref.ifT ) 
are shown by the symbols 
with the error bars, while 
the straight lines represent 
the calculated values. 





FIG. 9. Plot of the stopping powers (dE/dX) of and 
ions ys/E in UF 4 . The e 3 q)erimental points from 
Bridwell and Moak (Ref.lT-) are shown by the S 5 mibols ^ 
for and 5 ^®Br. The calculated values are shown 

by the straight lines ( for ^®Br and — for . 
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FIG. 10, Plot of the stopping powers (dE/dX) of 
Ions vs 'fS in C, Al, Ni, Ag, and Au. The symbols show 
the ejqjerimental values from Brown and Moak (Ref.^$) 
and the straight lines show the corresponding calculated 
values. The experimental errors (±10%) are not shown. 
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IV . 3 SIMPLER FORM OF STOPPING- POWER BQUATIOM 

A simpler^ but somewhat less precise, stopping power 

2 

equation can be obtained from Bohr' s classical equation for 
the energy-loss of a swiftly moving charged particle i.e., 

Eq . ( 3 . 6 ) : 


dx 


4 

mV^ 


„ ■ 1.123 mV^ 

nZ In 'Y- 


. . . (4.31) 


where to geometric mean value of the cyclic 


frequency of the electron in an atom of atomic nuirtoer 7^^ 
averaged over all the orbital electrons. Because of the non- 
participation of .many of the orbital electrons of a heavy 
medium in interactions with partially stripped heavy ions and 
because of a continuous increase in the number of such non- 
participating electrons with decreasing ion velocity, it is 
not meaningful to use Eq. (4 .31) with £5 in the logarithmic term 
Instead, one may retain the summation sign: 


dx 


4 TT e'^ ^n 

mV^ 


Z 

s 


In 


1.123 mV^ 


. . . (4.32) 


where the suiimatlon has to be carried out over all the s 
orbital electrons of the mediuo and (*1^ rs th.- cycl 

e the s'* orbital electron. Ihe logarithmic term 

qXiency of the s i o 

_eff „ .prz > , I =htd=-mir 

may be sirrplif ied by putting - ^ s s 2 s 



and = a /H , „here e and m ^ the electnenlc charge and 
™a3s, Is the ionization potential of the electron with 
orhital velocity and n is the nn„i.er of aton. of the ™edl» 
per unit volume. Poliowing Bohfs^ procedure, the sunuaation 
may be replaced by integration, with the upper limit u^= u 

where U^' makes the logarithmic term zero. Ohe logariLalo 
term in Eq. (4.32) is 


s U) 


= I In 5 2dl- ^23)V^ 

® lf(Z' ) u ^ 

J- S 


r fin i-ii^ldn(„ ) 

U =0 ^ ^ ^ \ TT ^ I s 


f Cz^) 


L 2.246. > 

''' 


. . . (4.33) 


Eq. (4.33) when used in Eq. (4.32) leads to the following 
energy loss equation: 


^ _ 8 (2.246) TTg^ 


n 


dx 


m V 


f(2,) 


{fCZj)]- 


3/2 V . 


(4.34) 


which on integration gives the following range-velocity 


rel ation 
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in which all the symbols are identical with those in Eq.(4.23) 
Since Eq.(4.23) predicts heavy ion ranges with accuracy, a 
comparison between Eqs. (4.23) and (4.35) shows that the 
latter would be useful if 



• • • ( 4r aSS) 


Actual computation of the above ratio shows that for all ions 
with 30{z^<45.5 it remains close to l.o5 while for Z^^45.5 
it varies from 1.055 to 1.07. Thus, adjustments of Eqs. 
(4.34) and (4.35) by a factor of 1.06 leads, to the following 
simple equations, which are about ±1% less precise than Eqs. 
(4.23) and (4.30) : 


i = b.4 

2 


R = 


809f(Z2) f 


...(4.37) 


, . . (4.38) 


Further Eqs. (4.23) and (4.30) may be somevj’hat simplified by 
noting the following; 


^ .7622 [fCz^)' 

j + f ( Z2)j>%r.20 + 0.5 Z^ 

for Z2^45.5 

. . . (4.39) 

|4.7622]f (Z^^)' 

j^/^+ f(Z^)|'=: 0.938. Z^ 

for Z^(45.5 

... (4 .40) 


with the introduction of ■~1% error. 
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IV, 4 DISCUSSION 



Equations (4.23) and (4.30) are based on the integra- 
tions, represented by Sq. (4.17) , performed with >.^1. For 
partially stripped ions this imposes the condition; 



for the applicability of Eqs. (4.23) and (4.30) . On the other- 
hand, Eqs. (4.37) and (4. 38), based oh purely classical treat- 
ment,^ require %^^1 or 2f (Z^)^^ 1 for their validity. Only 
in the case of U ion for which 2 f(Z^)=»9, this condition 
may be considered as fulfil'led. However, the calculated 
ranged of light fission products as also the stopping powers 
of "^^Br ion agree with the corresponding experimental data 

and 30 has been taken as the lower limit for the use of 

^9 14 

classical Eqs. (4.37) and (4.38) . For ions like Be and N, 
classical equations do not hold since for these ions the 
condition79) l is not fulfilled. For light ions C^Be ,^^N) 
the condition represented by Sq. (4 .41) is satisfied and 
Eqs. (4.23) and. (4.30) are valid. Ihus, Z^ = 3, where >- is 
just greater than unity, has been taken for the validity of 
Eqs. (4.23) and (4.30). 
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In the case of the media# the lower limit of are 

obtainable from Eqs.(4.7 ) and (4.33). If U ' is the value 

s 

of makes the logarithinic term zQro and if the 

orbital velocity of the K--shell electrons in the atom of the 
medium is Uj^, then for U^' ^ Eq, (4*.7 ) loses its signi- 
ficance* Use of this criterion requires that 

2 ...(4.42) 

This condition is obtained from the second integral of 
Eq. (4.7) # i.e.# 


U 


P , In } *1 an<U ) = 3 ^“1'^ 


zV^ V 


jd(U^) 


I . 


From the above equation the logarithmic term vanishes if 

U = and this leads to the condition given by Eq. 

s 

(4.42). Using 7.= 2f(Z^) in Bq.(4.42) one obtains 




... (4.43) 


as essential condition for the validity of Bqs. (4,23) and 
(4.30) . Similarly the classical treatment would be valid if 


U. 


K 


S ^ 2.246 7 

^tf(Z^‘ 


1/2 




...(4.44) 


V 
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This leads to the following condition regardirxg the ion 
velocity Vs 

v/u ^£51 1'/" 

^ K 1^2.246 J • <4.45) 


Eq.(4.44) is barely satisfied in the case of fission products 
in carbon as the medium. The fit of our range-velocity equa- 
tion in the case of Be may be either fortuitous or due to the 

relatively small contribution of the higher values of U to 

,s 

the total integral. In the case of hydrogen and helium the 
agreement appears to be purely fortuitous. In the case of 


heavy media, U, may greatly exceed the ion velocity at which 


the ion becomes completely stripped of its orbital electrons. 

* » 3IL sff 

Bohr s limitation that Z should always be less than 
— ^ for the applicability of Eqs. (4.23) and (4.30) does not 


seem to apply here. The assumption Z 


eff / ^ 


that 


./ 1 o 

^ 2 




requires 


. . (4.46) 


where V is the velocity of the ion, Z^ its nuclear charge 
and f(Z^) has been defined earlier. ; 

Eq. (4.46) is not always , fulfilled. There are velocity^^ 


regions in which the ion is partially stripped with 1 but 


z.v 

. 1 o 

Its velocity may exceed g'f (zT 
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IV. 4. 2 Deviation from Nidav' s approach 
12 

NidBY/ as a first approxiination, integrated 

Scj. ^4.11) t>etw6en the limits and on the assun^tion 
that at the moving ion undergoes large-angle scatterings 


and the resultant penetration in the beam direction is negli- 


gible. Niday observed slightly lower ranges in uranium by 

" using lead rather than aluminium to catch the fragments. To 

explain this discrepancy Niday assumed that significant 

scattering occurs along the path of the fragments in heavy 

media at V , the ion velocity at which the collision diame- 
c 

ter becomes equal to the screening radius and it was taken 

the 

as the lower velocity at which^nuclear collision becomes 
predominant. is obtained from the relation 




(Z^2/3, 2^2/3) 1/2 


..(4.47) 


where Z and Z are the nuclear charges of the ion and the 
1 2 

tli0 ■ ■ 

stopping atoms / is/^ reduced mass of the system and a^ is 

the ' first Bohr radius' of the hydrogen atom. Using as 
the lovjer limit for the velocity, one gets from Sq. (4.11) 


RCmg/cm'^) = 




127.3 f (Zj) [4.7622(kz/j^5 ^ 


1/3- 


. (4.48) 
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Calculations with as the lower limit, using the 
above equation, show that the values of k vary from 1.07 to 
1.12 as one goes from the heavy to the light media, (i.e., 
from uraiium to carbon and beryllium) . Using the new value 
of k instead of k=l in Eq. (4.5), it is found that in most 
of the cases (Al, Cu, Ag, Ni and Au) the agreement between 
the calculated and the experimental stopping powers tends to 
become worse but gives a better agreement between the calcu- 
lated and the experimental values of the stopping powers in 
carbon and beryllium. To optimize the agreement between the 
calculated and the experimental values of both ranges and 
stopping povjers, has been taken as the lov/er velocity limit 
for the range calculation. 

IV . 4 . 3 Stopping powers in carbon and beryllium 

The calculated stopping powers in carbon and beryllium 
have been consistently lower than the corresponding experi- 
mental values. On the other hand, the ranges of fission 
products in both carbon and beryllium, calculated with Eq. 
(4.23), are in good agreement with the corresponding experi- 
mental values. No attempt has been made to explain this 

■ • „22 ■ 

inconsistency but both Bridwell and Moak and Booth and 
Grant"^^ have mentioned that oxidation during thin film pre- 
paration by vacuum evaporation in the cases of both beryllium 
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and carbon and the tenacious retention of moisture in the 
case of Cctrbon may lead to unassignable errors in the 
measured stopping powers. 


iV,4.4 Range in the most favoured direction of r 
single C3rvstal 

If one assumes, that electronic collision is the only 
form of interaction by which the ion loses its energy till 
it stops, which may be approximately true in the "most favo- 
ured" direction in a single crystal, then the range-velocity 
equation, Eq. (4.23), should be modified to the following form 
in the case of a single crystal 


R(mg/cm^) = 


A, A V. 
1 , 


127.3 fiZ^) [4.7622j^f (Z^)|^/^ + 


...(4.49) 

H0nc0/ foir tli0 sBJTiG ion ox fission pxoduct/ bli 6 xstio of ttio 
ranges in polycrystal ^ mono-crystal ^n the 

most favoured direction) would be given by the ratio of Sqs. 
(4.23) and (4.49) . 


mono 


^ V. *1 
1 ' 

_ ^poiy J 

cal 

v,-v 

1 o 


. , . (4 .50) 


i 
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The validity of Eg. ( 4 . 49 ) can be ohec^^d bv c • 

^mono ^ comparing the rati, 

Vly Eg- (4.50) and from the e ■ 

^ ^ ® experimentally 

measured ranges in the :»st favoured direction of ■ 

ticrion of a single 

crystal R and in -hh^ ^ t 

mono “ ^ poly crystal R 

poly * 



CHAPTER V 

STOPPING- POWER EQUATIONS FOR 
COMPLETELY STRIPISD HEAVY IONS 

Bethe's energy loss equation is based on the assumption 
that , i.e. , 

V 

1 , ...(5,1) 

a condition which is always fulfilled in the case of low mass 
ions i.e./ protons and alpha particles , above 0,4 MeV/amu and 
Bethe’s equation j[Eq.(3.l2)J can be iised for stopping power 
calculations. In the case of heavy ions Eq,(3.l2) can be used 
if the velocity of the ion is high enough to satisfy Eq,(5.1) , 
but in practice this condition is not always fulfilled. Ihere 
are certain velocity regions in which » 1 and in this parti- 

4 

cular velocity domain neither Bethe’s equation nor Bohr' s 
classical equation^ is applic^ls* Eq. (4 .30) which was 
derived for the partially stripped ions with does not 

hold for completely stripped ions even though may be greater 


than unity. 
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It would be shown in this chpater that by a careful 
choice of the limits of integrations appearing in Eq. (4*1)^ 
it is possible to calculate the stopping powers of heavy ions/ 
including proton, at an energy which is high enough for elec- 
tronic collisions to be the precominant mode of energy loss> 

The derivation, calculation, and discussion of these are 
presented in the following sections. 

V.l DERIVATION 

V.1.1 Energy-loss equations 

The logarithmic terms in Bohr' s energj'-loss equation, 

Eq. (4.1), may be written in the form of integrals as discussed 
in Chapter TV. 

s 

• (5,2) 

■ V. ■ ■ 

where ^ = ~ and5C-= 2z .^. Ihe second integral in the above 
® s 

equation may be evaluated in two parts, as shown in Sec. IV. 1.1 

. . . 

of Chapter IV, corresponding to the limits within which 
is less than and greater than unity respectively. Indicating 

the three integrals by J 2 , J 2 a one gets 

J = *^2 ''' '^3 ■ ' 


( 5 . 3 ) 
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U =U ' 
S/ s 


where J, = 


.L -IVM-} 


c3n(U ) 
s 


. . . (5.4) 


U3=2V-X.- 


J. = 


7 


In 


U =0 
s 




» » • (5.5) 


and 


J. = 


U =U " 
s f s «• 

-1 r ^(u„) 

Ug=2V>'^ L ® « J ® 


— V. 


. • • ( 5 . 6) 


Each integral raguires careful examination for fixing the 
upper limit of integration corresponding to different actual 
physical situations* TWo cases may be considered depending 
on the values of ^ . 

(A) *X>1 

(a) 

n 2 V 

— / the first integral may be written 
. s 
as 

= J ln((-P) dn(U^) , ...(5.7) 

all U L s J 

S ' 

where the upper limit of the integration is given by that value 
of for which the logarithmic’ term vanishes. Ihis puts the 

condition 2 vNu % for the logarithmic term to be positive. 
Further, since the maximum value which U can assume corresponds 
to the orbital velocity of the K-shell electrons of the me dium^ 

the logarithmic term would be always positive if 
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Z^V 

V 2 




. . . (5.8) 


where the velocity of the K-shell electrons is assnmed to 
be ZjV^. Putting i^= j mu^Z, In Eq.(5.7), Jj thus mey be 
written as 


j 


= J In ) 

all U ■ I 7C ® 

s s 


= in 


2mV" 


. . . (5.9) 


■Ml 

where I is the ionization potential of the s orbital 

O'. 

electron of the medium of atomic number and I stands 
for the geometric mean ionization potential of the medium, 
obtained by averaging over all the orbital electrons. 


Z^ In I = J In I dn(U ) 

2 all U s s 

S 


...(5.10) 


As far as the integral J 2 is concerned, the logarithmic term 

lln(~)^^ would be positive for all values of upto Z^V^, 

^ V® ZoVq V ^2’^ov 

if V / — — . Since, the condition — -a- has already 

^ \ ^2^ 

been assumed, V/ — ^ 'is automatically satisfied for'Xyi. 
Using I' = ^ mU ^ and Eg. ( 5 . lo) , .in Eq. (5.5), J, may be 
given by 


all U 


( ihj(f>^<ah(v 


L s 


= Z 2 in 


2mV" 


...(5.11) 
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Further, since the upper limit for u in J is 2V>r^, the 

s 2 ^ 

third integral becomes superfluous for Z V , 

^ 2 o 

since the lower limit for the integral already corresponds 

to the maximum possible value for U , Thus, the combination 

s 

of and J^, given by Eqs.(5.9) and (5,ll),and substitution 
in sq. (4.1) lead to the following expression for the energy 
loss: 

r2 


dE 4irz^e^ „ „ , 2mV 

S - " ^2 1° 


/ 


mV" " Ipt 

where all the symbols have been defined earlier. 


...(5.12) 


(b) 




Z_V . Z^V z^v 

since y^l, y^~ 2 ~ ' and and 


are evaluated as follows: _ ■ 

Jt requires an upper cut-off value U =U ' , since the 
J. s s 

logarithmic term ^In (§^) would otherwise assume negative 

values for some values of U_ / Z«V_, Hence from i:;qs.(5.4) 

s / o 

dtJ 

and dn(U^) = f(Z 2 )y:^/ 


U ’ 
'S 


= r 

u =oL 


o 


2% 

2f (Z-,) 

Ife) = 

2 

S J • ^ 

"O , 


4f (Z2)V 


V 7^ 

o 


U V 

r\ ' 


u 


2V/(j 


dU 


. . . (5.13) 
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Because of the condition that the second Integral 

Jj, given by Eq. (5.5) does not- include all possible values 
of u upto Z V ; 

S 2 O 




2f(Z^) 

2VX"^ 

V 

' 1 ln( 

0 

0 

4f(Z,)V 

~ ^"o7^ 

(l+ln%) , . 


2V^ 


The third integral is given by 


J., = 


U^=2V;7C^ ^ ® J 


3f(z2) ys 

V 


% 

2V X ^ ® J 


. . (5.15) 


-1/'^ 

The above integral vanishes for U = 2V'X/ . Since the 

maximum value of U is assumed to be Z^V , if Z„V exceeds 

. s 2 o 2 o 

-1/3 

2V pC then an upper cut-off value has to be used for evalu- 
ating the integral because for 2V '^ 2^0 is physi- 

c ally meaningless to evaluate the : integral in the domain 

2V V,U > Z V . The upper cut-off value for the third 

^ sV 2 o 

integral may be obtained from a consideration of the number 

n (U ) of orbital electrons with u\ 2V X which interact with 

' s ■ , ■ -S'- 
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the ion at a given velocity v. m the actual evaluation of 
J3 It is convenient to consider the contribution of the t™ 
K-shell electrons separately from that of the other electrons 
with \ This procedure avoids the complications 

arising out of the relativistic velocities of the K-shell 

electrons in heavy media aa ttoI i +.v ■ 

avy Hieaia as well as the inaccuracy involved 

in using the expression for n(up, given by Eq. (4.16), in 

the case of the inner-most electrons. Eq.(4.16) for n(u ) 

suffers from the drawback that it does not lead to the coLct 

velocities of either the innermost or the outermost electron. 

Further, the range and the stopping power data, on which the 

validity of Eq. ( 4 . 16) is based involve only a few outer 

electrons in the case of heavy media and may not hold for 

the deeper-seated electrons. Taking these limitations into 

considera cion, may be written as 


(Z -2)V^ 

U = 2 _p 

s ,.f(Zj 




Ug=2VX^ 


.dn(U^)' + 2 In j(f~) X. 


-1 


2 o 


(Z--2)V^ 
2 o 


3f(Z2)V 


V 


I ; H x- j 


dn(U ) + 2 In 
s 


^ j 


In ^ (z - 2 ) 


2Vf(Z ) ,1/3 

> + 3 ( 23 - 2 ) + 




•^2V >r^(ln )C^^^-l)j 


+ 6 in - 2 InTt] . . . (5.16) 


22''o J 
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J is given by the sum of Eqs. (5.13) , (5.14) and (5.16) which 
when incorporated into Eg. (4.1), gives 


dx 


2Trz^e^n f, 2Vf(zJ 

[3(Zj-2) + 3 (Z 2 - 2 ) in + 6 In 5 ^ 


2 o 


2f(Z^)V 

+ _____ _ 2 ^ ln> 


-J / ...(5. 


17) 


V 


all 


where 9 C= 22 ~ and^the other symbols have been defined 
earlier. 


(c) v(, 


Z V • I/O 
2 o 


If v( 


Z2V ^ , 

2 5 *^ then both and would require 


appropriate cut-off velocities U since all the orbital 

■ S'. 

electrons can not have interactions with the rroving ion. 
This case has already been considered in Sec. TV. 1.1 of 
Chapter IV. The energy loss is given by 


dx 


2 4 f(Z_)V 
4^2 e ^ 2 


mV 


V 


[ 


3 


] - 


(5.18) 


(B) For%(^ 1 

For %X l/ since is u^^ in the integral J^, both 

the integrals and becorne identical . Thus, J is given by 
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^ 2 , 1 in V dn (U) 

all U s 

s 


. • .(5.19) 


The third tategral does not come Into the picture, because 
for >-^1, (-jj) IS always less than unity, as shown In 

Sec.III.2.5 of Chapter III. Prom Eg. ( 5 . 19 ) 


2 Jin (^)^ ,dn(U^) 


= 21^ in 


* . . ( 5.20) 


where I is given by Eg. (5.1o) . Thus, the’ total energy loss 


becomes. 


_ dB _ e n „ , 2mV^ ■ . 

dx “ 2 ^2 ...(5.21) 

This eguation is identical with Bethe' s eguation except that 
I stands for the mean ionization potential as given in Eg. 

(5.21), while that in Bethe* s eguation I stands for the mean 

3 ' ' 

excitation potential. Bohr showed that the mean excitation 
potential can be put egual to the mean ionization potential 
of the electrons with sufficient accuracy. 


V.1.2 


!s equations for completely stripped heavy ions 


Except the proton and the alj*ia particle which have 
un(iex all 0031(111:10313 as i^^ntioned earlior/ all othex 
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heavy ions would have » 1 l„ ,,, inoonpletely stripped 
state, ™ay have X> 1 certain velocities in the completely 

3Xid a.t a.DDiroi')ri“'t'i=> Vinrrn t i • 

dppropricite hign velocities would have 
Since completely stripped ions are considered, the 
condition which should be fulfilled is that at the minimum 
velocity, the ion should have its charge number 2 
equal to the nuclear charge Z^. Hence, 


eff . V - 

Z = f(Zj) - = Zj , 


. . (5.22) 


where f(z^) is given by Eqs . (4 . 24) and(4 .25) and V^= e^fi. 
This ^ gives 


Z,V 


- V 

_ij2 


f(Zj) ' 


. . (5.23) 


as the essential condition for the ion being completely 
stripped. For such ions the energy loss equations are given 
by Eqs. ( 5 . 12 ), (5.17), (5.18) and (5.21) with 2 = After 
substitution of the values of the physical constants in the 
energy loss equations, one obtains the following simplified 
relations: 

(A) 7.^ or V ^ 2 Z^V^ 

ZoV. 




- ^ (JJg31S!s!) = 63.65 log - ...( 5 . 24 ) 

cbc ma t 


mg 
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2 o 


TCN \^2''o V3 Z,V4 3/4^ 


^ (MeV__cmf. 
dx mg . 


13.79 — 3(Z -2) + 3(Z„-2) In 

V L' ^ ^ 


zvfCz^) 

'iZ.- 2 )V^ 
I O 


2y 2f(z ) V 

+ 6 In - - ' + —2 — 

2 o o / 


^2 Injt; 


. . . (5.25) 


/^2^o 3/3 ,Z 

(c)V<— 5^ orv(-i_| s 


f(z^) Z 2 


dE (MeV_cm:_) _ .2 68 ^ r 3 - - 

dx ^ mg ^ - ^2.68 ^ |^3 7 ^ + % 


'J ■•<5-: 


(B) For 1 or v)2Z V 


^ ( Mey cmf ., ^ (ii^Sy!) , 

dx mg a v2 - I 


..(5.27) 


where I is the mean ionization potential in -units of eV, V is 

the velocity of the ion and is the electron velocity in the 

first Bohr orbit, both in units of lo® cn/sec; £(^2) given 

^o 

by Eqs .(4 . 24):^d (4 .25) and 'X = 2Zj — . 
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V-1.3 Range- enercfV relationship 

The total range may be obtained from the energy loss 
equation by simple integration. Two cases arise depending 
on the values of'^. 


( A ) '^±1 

For 1 there are three energy- loss equations/ which 
are valid only in certain energy domains. Using S = V^, 


V 


dE 


M^VdV and,?C= 2z — , in Eq. (5.12)/ one gets 


M- V^dV 

dR = dx = (-^) — ^ -pr- — 

47re^ nz^ Z In ( ^ - ) 
2 2 V" 1 


. . . (5.28) 


where M^= A^in^, m^ being the nucleonic mass. Tliis equation 
is not valid at velocities given by mV^^V^I, since, the loga- 
rithmic term in Eq. (5.12) becomes negative and loses its 
significance. Therefore Eq. (5.28) can not be integrated 
over entire velocity range. But it can be used to describe 
the difference in range between any two hi^ veloci- 

ties and V^. Eq. (5.28) may be written as 


dR = (- 


m 


A^m V'" dV 
1 o 


47re' 


. . . (5.29) 


3 

"2 2 - 


I 

iW 


nz^ Z„ - In I ( _ V 

z V I 
o 


s 2/3 ,,2 


] 


1 ^' 
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Siibstituting 


[ 


(-^) 
zV I 
o 


2/3 


V 




• « » ( S * 30} 


in Eq.(5.29) and integrating between the velocities and 
V^/ one obtains 


R — 

(2 _ A-m z V I 4/3 /'2 

) dR - ( 4> — i taC •••<S.31) 


or (R^-Rj) 


„ V 4/3 AtIH 
/ Xo 

^ 2/3. 


F4/3 


12'rre 


n z " 


j^Ei(lnC2)-Ei(lnC^)j, 


. . . (5.32). 

where Ei stands for the exponential integral and or C 

is given by Eq. (5.30) . Eq. (5.32) holds for all velocities 
. z V I 1/3^* 

(B) ><■ 1 

Eq. (5.21) on integration gives 


(Rj-Rl) = ( 


m 


- 2 A-V^ p 

) (-) — i Eidn d,)-Ei(ln d ) 

m 2_ i J- • ^ 


4" 'm' ^ 2- t 

32 Tre ^2 


. . (5.33) 


2 2 


where d,/ d. or d is given by d = (-^^) / Eq. (5.33) as valid 

■ 2 ; ; ■ . - x/2 ^ ■ 

for all velocities * 
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Thus from Eqs. (5.32) and (5.33) the range difference 
at two different velocities can be obtained, using the expo- 
nential integrals. 


relat i onship for completely striDoed 
heavy ions — 

Substituting z = and the values of the constants 
in Eqs. (5.32) and (5,33)/ the range-energy relationship 
becomes^ 


(A) For 1 

(mg/cm^) = (3 .494xlo'*^) 

f 


A I 

V_„ 


2^2/3 


jEidn C^) ~Ei(^ln C^)| ... (5.34) 


where A^, A^ and Z^, Z^ are the mass numbers and atomic numbers 


of the ion and the stopping medium respectively/ I is the mean 
ionization potential of the mediiim in eV and C^, or C is 
given by 

1 / n.^q -s'" ^ 
c 


^ 11^39 
2V_I Zt 


...(5.35) 


8 


In Eq. (5.35) V is the veioeity of the ion, and V^=2. 18x10 cn/sec. 


.8 


both in units of iO cirj/sec • 
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(B) For^C l 


5 Aj AjI2 


[Rj-hJ (.„a/cn.2, = (1.453X10 -) — -f-jadn d,)-a(i„aj 

Z- 2 - ■ I_ 


• . . (5 . 36) 


where d^, or d is given by 



...(5.37) 


V.1.5 Mean ionization potential 
.51 

Mukherji has derived an expression for the calcula- 
tion of the mean ionization potentials of elements as follows. 

The geometric mean ionization potential I of an atom 
of atomic number Z is given by 

Z In I = I In I ... (5.38) 

s s 

where is the ionization potential of the s^ orbital 
electron of the atom and the summation includes all the Z 
electrons. Adopting Bohr's^ procedure, the sximmation can be 
replaced by an integration. 

I in I = j |ln mU j dn(U ) ...(5.39) 

s s all U 1 2 ® J 

S: / 

where m is the electronic mass, being the velocity of the 
s^^ orbital electron. The upper and ttie lower limits of 
integration should be the Velocities of the innermost and the 
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outennost electrons in the atom respectively. Since the 
expressions for n(U^) , given by Eg. (4.16) do not lead to 
correct values for the velocities of either the innermost dr 


outermost electron, in evaluating Eg. (5.39), it is convenient 

separately 

to consider the contribution of two K-shell electrdhs/from 


that of other (Z-2) electrons. 

U =U ' 
s r s 

J ln(^mu3dn(U) = 
U =0 ^ S s 

s 


(Z-2)V 

O 

f (z) 

J ln(^ mU dn(U ) 
z s s 

o 

+ 2 In (13.6) Z^ 

...(5.40) 


The lower limit of integration U =0 is an approximation made 

s 

' ■ 3 ■ 

earlier by Bohr. In the case of lii^t elements, the ioniza- 
tion potential of each K-shell electron has been taken 
2 

as 13.6 Z eV and this seems to hold in t±ie cases of very heavy 
elements like uranium, as the electron binding energy data of 
Hagstrom et al.'^ show- With dn{U^) = dU^ f3X)m Eq«(4^16)jr 

O' 

Eg. (5.40) gives ■ . 

Z In i = (Z-2) ln|^^|^|^^j^.l3.6 + 2 In 13.6z2, . . (5 .41) 

where I is in electron volts, Z is atomic number of the atom 
and f(Z) is given by Egs. (4.24) and (4.25). 
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V.2 CALCULATIONS AND COMPARISON WITH EXPERIMENTAL VALUES 

Limited number of experimental are availa- 

ble on the energy loss measurement of heavy ions in high 
energy regions, in which the ion is cornpletely stripped. Ihe 
condition for the ion being completely stripped is given by 
Eq. (5.23) . Table V.l shows the minimum energy Em for a few 
ions at which the ions become completely stripped. Ihe last 
column of the Table V.l gives the minimum energy of the ion 
required for %^1. Once the requirement for the ion being 
completely stripped is fxalfilled, the velocity of the ion 

must also meet either of the conditions V /,YZ-Z’ V and 

Zl/4z3/4v <^'120 




for the applicability of Eqs. (5.24) and 


(5.25) . Table V.2 gives the values of the minimum energies 
E. and Bn required to fxilfill the conditions VvY2-Z.V and 

^ Z,V222 374 V 12 0 

Vj;j ^ for different incident- ion stopping-medium 

J2 

combinations. 


V.2.1 Stopping-powers 

Hower and Pairhall'^^ have measured the energy degrada- 
tion of ^Be ions in gold and aluminixim foils. Table V. 3 shows 
the e 5 <p e rime nt ally obtained values of the stopping powers and 
the corresponding calculated values with Eq. (5.26) « Experi- 
mental stopping powers at different energies were obtained 
in the following manner. Prom the incident and eitergent 
energies E^ and E^ of the ion before and after passing throu^ 
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TABLE V.2 


The minimxira projectile energies S- and 



required to fulfill the conditions 




2 o 



All the energies are in xjnits of MeV/amu. 


Ion 



S 

TOP 

P I N 

G M 

EDI 

U M 





^2 

"" ^2 

El • 


„ Au 
^1 

E j 

2 .. 

Be 

0.56 

0.79 

1.16 

1.29 

1.30 

1.38 

3.67 

2.77 

17.39 

7.8: 

B 

0-62 

0.99 

1.30 

1.61 

1.45 

1.73' 

4.10 

3.46 

19.44 

9,7^ 

C 

0.68 

1.19 

1.42 

1.93 

1.59 

-2 .03 

4.49 

4,16 

21,30 

11.74 

N 

0.74 

1.38 

1.53 

2.25 

1.71 

2,42 

4.85 

4,85 

23.01 

1 

13. 7d 

0 

0.79 

1.58 

1.64 

2.57 

1.83 

2.77 

5.19 

5.54 

24.60 

15.65 

F . 

0.84 

1,78 

1.74 

2.89 

1.94 

3.12 

5.50 

6.24 

26.09 

17.61 

Ne 

0.88 

1,98 

1.88 

3,22 

2.05 

3.46 

5.80 

6.93 

27.50 

19.57! 
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TABLE V.3 


Energy (MeV) 

(dE/dx)Exp?^ 
(MeV cm^/mg) 

(dE/dx)Cal. 
(MeV cm^/mg) 

9.75 

1.13 

1.11 

13.85 

0.90 • 

1.02 


a foil of thickness t, the experimental stopping power is 
(E.-E ) 

obtained as — . The corresponding calculated valiaes 
have been obtained from Eg. (5.26) using a value of V corres- 
ponding to the mean energies experi- 

mental data on dEJ/dx could be obtained in the case of gold. 


while in the case of aluminiiam no reliable data could l^e 
obtained from the energy vs thickness plot. In a similar 
manner the energy-degradation data of Northcliffe have been 


treated to obtain the experimental and calculated stopping 
powers for the and ^°Ne ions. Tables V.4 - v,7 

show these results in the case of aluminium. The mean ioniza- 


tion potential X for alumini^ 2 m has been calculated with Eg. 
(5.41) and is 151.4 eV. Kelley et al. have measured the 
energy loss of and ions at energies of 2-10 MeV/ 

nucleon in 94.6 thick silicon detector. The experimentally 
obtained values of the stopping powers along with the calcula- 
ted values, using Eqs. (5.24) and ( 5 . 27 ) , are shown in Tables 
V.8 -Y.IO. For silicon 1= 158.8 eV was used m ^ 
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TABLE V.4 

Experimentally obtained values of the stopping-powers of 

14 . 97 

N ion in aliminium along with the calculated values 


^ ,MeV. 

Energy 

( dE/ dx) Exp . 
(MeV cm^/mg) 

( dE/dx) Cal . 
(MeV cm^/ mg) 

Equation used 

9.53 

1.69 

1.76 

Eq.C5.27) 

8.30 

1.81 

2.03 

Eq.(5.27) 

7.46 

2.23 

2.13 

Eq. (5.27) 

6.34 

2.46 

2.42 

Eq.(5.27) 

4.87 

2.95 

2.97 

Eq,(5.24) 

4.12 

3.22 

3,31 

Eq.(5.24) 

3.33 

3.62 

3.76 

Eq.(5.24) 

2.66 

4.42 

4.28 

Eq.(5.24) 

1.85 

4.87 

5.03 

Eq. (5.25) 

1.41 

5.08 

5.11 

Eq. (5.26) 
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TABLE V.5 

Experimentally obtained values of the stopping-powers, of 
^^0 ion in alumini'um along with the calculated values^*^ 


„ /MeVv 

Energy 

(dE/dx) Exp. 

2 

(I'leV cm /mg) 

( dE/dx) Cal 

2 

( MeV cm /mg) 

Equation used 

10.25 

2.08. 

2.14 

Eq. (5.27) 

9.11 

2.40 

2,38 

Eq. (5.27) 

8.24 

2.64 ■ 

2,57 

Eq, (5,27) 

7.47 

2.92 

2.78 

Eq. (5.27) 

6.01 

' 3.11 

3.28 

Eq. (5.24) 

4.80 

3.73 

3.79 

Eq. (5,24) 

3.65 

4.31 

4.49 

Eq. (5.24) 

3.04 

4.8 3 

4.99 

Eq. (5.24) 

2.30 

5.77 

5.86 

Eq. (5.25) 
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TA3LB V.6 

Experimentally obtained values of the stopping powers of 
ion in alxminixom along witli the calculated values^"^ 


Energy ( J 

(dE/dx) Exp, 
(MeV cm^/mg) 

(dE/dx) Cal . 
(MeV cm^/mg) 

Equation used 

9.60 

2,85 

2.89 

Eq.(5.27) 

9.04 

2.98 

3.08 

Eq.(5.27) 

7.92 

3.14 

3.35 

Eq.(5.24) 

7. 40 " 

3'. 49 

3.51 

Eq. (5.24) 

6.47 

3.82 

3 .83 

Eq.(5.24) 

5.49 

4.53 

4.31 

Eq. (5.24) 

4.56 

4.73 

4.80 

Eq.(5.24) 

3.57 

5.94 

5.68 

Eq.(5.24) 

2 .48 

6.60 

6.82 

Eq.(5.25) 
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TABLE V. 7 

Experimentally obtained values of the stopping-powers of 
20Ne ions in aluminium along with the calculated values27 


r, /MeVq 

Energy 

(dE/dx) Exp. 
(MeV cm^/mg) 

(dE/dx) Cal. 
(MeV cm^/ma) 

Equation used 

9.92 

3.51 

3.46 

Eq.(5.24) 

9.23 

3.86 

3.65 

fl 

8.13 

4.00 

3.98 

11 

7.19 

4.21 

4.32 

« 

6 . 20 

4.86 

5.02 

11 

5.38 

5.04 

5.22 

» 

3.50 

6.51 

6.74 

11 

TABLE V.8 


Experimentally obtained values of the stopping-powers of 
ion in 94.6 /^m thick Si-detector along with the cal- 
culated values using Eq. (5,27)2® 


E. 

1 

(MeV) 

(MeV) 

— As 

S=E.- ~ 

(MeV) 

(dE/dx) Exp. 
(Kev cmVmg) 

(dE/dx) Cal. 
(KeV cm^/mg) 

120.0 

32.0 

lo4 .O 

1451 

1432 

109.2 

35.0 

91,7 

1588 

1580 

102.0 

36.0 

84,0 

1633 

1691 

92.0 

38.0 

73.0 

1724 

1888 

84.0 

49.0 

60.5 

2223 

2179 

78.0 

54.0 

51.0 

2449 

2482 
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TABLE V.9 


Experimentally obtained values of the stopping-powers of 

ion in 94 »6 ^^m thick Si— detector along with the cal— 
ci-ilated values^® 


E. 

(MeV) 

AE 

(MeV) 

f 

(MeV) 

( dE/ dx) Exp . 
(KeV cm^/mg) 

(dE/dx)Cal. 
(KeV cr^/mq) 

140.0 

42.5 

118.8 

1928 

1999 

127.0 

48.0 

103.0 

2178 

2237 

120.0 

50.0 

95.0 

2276 

2382 

110.0 

56.0 

82.0 

2541 

2673 

loo.o 

65 .0 

67.5 

2951 

3099 


TABLE V.lO 


Experimentally obtained values of the stopping-powers of 
^^0 ion in 94,6 Um thick Si- detector along with the cal- 
culated values^^ 


E . 

1 

(MeV) 

AE 

(MeV) 

v_E - ^ 
i 2 

(MeV) 

(dE/dx)Exp. 
(KeV cm^/mg) 

(d^/dx)Cal. 
(KeV cm^/mg) 

Equation 

used 

155.0 . 

56 .0 

126.0 

2630 

2736 

Eq.(5.27) 

140.0 

68.0 

106 .0 

3083 

3221 

Eq. (5 .27) 

132.0 

77.0 

94.5 

3491 

3436 

Eq. (5.24) 

122.0 ■ 

86.0 

79.0 

3899 

3856 

Eq. (5 . 24) 

ilo.o 

100.0 

60.0 

4539 

4564 

Eq. (5.24) 
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and (5.27), for the calculation of stopping powers. The 
value of I was obtained using Eq. (5.41) . 

V . 2 . 2 Range- differences 

Table V. 11 shows the range difference '.R^-R^) at two 

different energies (E^-E^) for and 

20 

Ne ions in aluminium. The ej<perimental values are obtained 
from the energy-loss data of Northcliffe , v/hile the calcu- 
lated values are obtained using either Eq. (5.34) or Eq.(3,36) 
depending on the values of % greater or less than mlty. Ihe 

exponential integrals were evaluated using the "Tables of 

50 

sine, cosine and exponential integrals". The uncertainty 
in the experimental ranges are not shown in the table. 

V.3 DISCUSSION 

The stopping power equations given by Eqs, (5.24), 

(5. -25) and (5.26) can be used for any completely stripped 
ion with %^1, while Eq. (5.27) is v?lid only for those ions 
for which X^l. Ihe condition for 7<.<: 1 for ions heavier tiien 
protons and alpha particles always fulfill the condition 
given by Eq . (5.23), i . e . , f or the ion to be in a compl etely 
g "tripped state. Ihe agreement between the calculated and the 
experimentally obtained values of the stopping powers is 
usually good. 
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TABLE V.ll 

Range-Differences for heavy ions in 


27 

altaniniiom in terras of (R 2 )-('R 2 ^ * 

Ion 

Energy interval 
(MeV/arau) 

Range-Exp . 

2 

(mg/ cm ) 

Range- Cal, 
(Mg/cm^) 

, Equation 
used 

“b 

9.73-2.92 

56.97’ 

54.36 

Eq. (5.36) 

“b 

9.88-2.48 

66.49 

64.03 

- Eq. (5.36) 

llg 

2.48-1.68 

3.53 

3.21 

Eq. (5.34) 


9.00-4.74 

36.72 

35.82 

Eq. (5 .36) 


4.74-1.65 

48.70 

46,61 

Eq. (5.34) 

1^0 

9.41-6.34 

25.06 

24.15 

Eq.(5.36) 

O , 

6.34-2.55 

41,49 

39,21 

Eq.(5.34) 


9.69-8.04 

14 .14 

13.51 

Eq. (5.36) 

19p 

8.04-3.00 

37.92 

35.98 

Eq. (5.34) 

^“^Ne 

9.35-3.17 

31.89 

30.52 

Eq.(5.34) 



chapter VI 


EFFECTIVE CHARGE 

The charge of a fast ion moving through matter fluctuates! 
as a result of the competition between electron loss and electron 
capture in collisions with the atoms of the stripping target. 
After the ions have made a sufficient nimber of collisions with 
the target atoms an equilibrium distribution of charges is estab-: 
lished which is dependent only on the velocity of the ions. Ibis 
equilibriiam distribution of charges does not change when the 
target thickness is. further increased. The mean equilibrium 
charge q (equilibrium charge) is obtained fTOm the equilibrium 
charge distributions. For the purpose of stopping power calcu- ^ 
lations, q is usually identified with the "effective charge” 

(Z ) of the ion which hashbeen described earlier. 

^ T^^ excellent agreement of the calcxiLated values of the 

stopping powers of various ions ( Be to U) with the corres- 
ponding experimental values in different solid media except 
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carbon and beryllium, seems to provide very strong, although 
indirect, evidence in favour of the assunption that 

= f(Z) ^ . ...{6.1) ^ 

O ■ 

where f(Z) is given by Eqs. (4.24) and (4.25). The fact that the 
light gases, as shown in Fig. 1, also obey Eq. (4.20), possibly 
indicates that Eq, (6,1) is independent of the physical states 
of the medium. In the following section a comparison of the 
measured equilibrium charges q and the corresponding values of 
calculated from Eq. (6.1) has been made to establish the 
validity of Eq. (6.1) . 

VI , 1 Comparison of the directly measured equilibrium charges i 
and the calculated effective charges 

Tables VI, 1 - VI. 6 show some experimental values of q 

9 

5 5— 57 

and ■g from some authors along with our calculated values 

O' 

using Eq. (4,22), where q and q are the equilibrium charges of; 

g s 

an ion after it emerges from a gaseous and a solid medium. In 
some cases the standard deviations associated with the measurec 
equilibrium charges are also listed in the t^les. The experi— 
mental set up for the measurement of q consists of the following: 

(a) A magnetic analyzer to select accelerated ions of 

xinique charge states. . [ 

(b) A stripper chamber which contained a foil and a yj 

gas cell. 
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(c) An snalyzsr to sspansits tho bsam ©msirg'ing' fjrom th© 
foil or the gas stripper into the various charge states. 

(d) A detection system to determine the relative intensi 
ties of the charge state conponents. In this way post stripper 
charge distributions and their corresponding pre-stripper charg 
can be obtained for different energies simultaneously. 

From the Tables VI. 1 - VI. 6 it is evident that for all 

— jF ' 

ions the agreement between q and the calculated Z is good ; 
ef f — 

but Z differs considerably from q^. This difference in the 
values of Z and g can be explained if it is ass^med that 

inside a solid medium the equilibrium charge does not differ 
appreciably from the equilibrium charge in a gaseous medium. 

' 58 ,"' . . 

BetZ/ in his review work, has agreed with the conclusion of 

Betz and Grodzins that inside a solid medium the equilibrium ; 

charge is possibly the same as that in a gaseous medium. The 

gas-solid difference has been attributed to the fact that a suff 

cient time elapses between inelastic collisions in a gas for ion: 

to de-excite to their ground states whereas in a solid many coll: 

sions take place while the ions are in the excited state from 

which electrons are readily stripped. The higher value of the 

equilibrium charge when the ion emerges from, a solid stripper ma] 

5o ' ' 

be attributed to the surface effects. 
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TABLE VI. 9 

Directly measured equilibrium charges (q^) in gaseous media 

for various ions along with the calculaued values of the effec- 

. ©f f 

-“tive charge Z at different energies 


Ions 

Velocity x lo"® 
crrv^sec 

He 

“2 

Ar 

Kr 

^eff 

(Calcd) 

“o" 

6,29 

- 

- 

3 * 5+0 #766 

- 

3 . 23 


■ 8.55 

- 

- 

4.8+0.766 

- 

4.39 


10,21 

- 


5.2+0.766 

- 

5,24 

20,, 

Ne 

2 . 65^ 

1.6 

1.7 

1.7 

1.5 

1.58 


4 .04^ 

2-2 

2.2 

2.5 

2.2 

2.41 


5.45^ 

- 

- 

3 » 3+0,766 

- 

3.25 


5.52^ 

2.7 

3.5 

3.5 

3*3 

3.29 


6 . 3 3^ 

3.4 

4.0 

- - 

- 

3.77 


8.25^ 

- 

- 

4.9+0.766 

- 

4.92 


9.94^ 

- 

- 

6.0+0.766 

- 

5.93 


a. Hubbard and Lauer, Phys. Rev., 98, 1814 (1955); 

b, Nikolaev et al., JSTP 12, 627 (1961). 


TABLE VI. 3 

23 31 40 

Directly measured equilibrium charges of Na, P and Ar 
ions in ^^(qg.) along with the calcxilated values of Z . 


Ions 

Velocity x 10 ® 
cn/sec i 


Z®^^(Galcd) 

23 ,, ^ 

Na 

4.5 

3.1 


3 lpa 

4.05 

3.0 

3.16 ■ 

40 ^ a 

Ar 

4.05 

3.2 

3.57 


a,' Nikolaev et al-, JETP, 12, 627 (1961), 
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TABLE VI. 5 

• 3 . *^9 

Directly measured equilibrium charges of S ions in 

air (q ) and carbon (q ) along with the calculated values 

efF ® 

of Z at different energies 


Energy (MeV) 

qg(air) 

^■f -f 

q (carbon) Z 

s 

(calculated) 

3.70 

3.45 

5.85 

3.84 

6 ,40 

4.65 

6.85 ‘ 

5.05 

8.40 

5.16 

7.45 

5.78 

10,00 

5.60 

7.80 

6.31 

12.70 

6 #50 

8.45 

7.11 

17.00 

7.50 

9.09 

8.23 

20.00 

8.06 

9.45 

8.92 

22.70 

8 . 50 

9.71 

9.51 

23.80 

8.65 

9.87 

9.73 

26.40 

8.90 

10.14 

10.26 

■27.20 

8.80 

10.25 

10.41 

29.90 

9.28 

10.46 

10.91 

32 . 60 

9.58 

10.66 

11.40 

a, Betz et al.. 

Phys, Letters 

/ 22, 643 (1966) , 
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VI. 2 Discussion 

A ^ 8,48,49,60,61 ,, ■ ^ . 

According to many authors . , the stopping 

power of a heavy ion is related to the stopping power of proton 

at the same Velocity and in the same stopping medium hy 


Q'p "p 2 

(Z (dE/dx)j^^ 

^ (dE/dx)p 


- ( 6 . 2 ) 


where the subscripts HI and p refer to the heavy ion and the 

proton respectively. Eq. (6.2) is obtained from Bethe' s equation, 

in which the logarithmic part is independent of the charge of the 

eff 

ion. Eq. (6.2) has been utilized to make an estimate of Z 
from the measured stopping powers and for obtaining the energy- 
loss and ranges of heavy ions through interpolations and extra- 
polations. Bohr pointed out that Eg. (3.12) is applicable if 

or 2f(Z^)^l, a condition which is not always firlfilled in 
the case of partially stripped heavy ions. Cn the other hand, 
if 2f(Z^)^l, which is marginally valid for most heavy ions, 

Eq. (4.37), which is derived from Eg. (4.32), can be utilized for 
stopping power calculations. Since Eq. (4.37) predicts stopping 
powers that are in fair agreement with the experimental values, 
a comparison between Eqs. (4.32) and (3.12) shows that the assunp- 
tion of the logarithmic term being simply velocity dependent may not 

.c .c 

be strictly valid. The values of Z obtained with Eq. (6,2) may, 
therefore be somewhat different from the actual values inside a 


solid medium. 



CHAPTER VII 


SUMMARY 

Using Bohr's semi-quantum mechanical expression for the 
energy loss of a swiftly moving charged particle a set of stopp- 
ing power and range-velocity equations have been dedxx:ed. These 
equations predict stopping powers with an average deviation of 
less than 5 % from the experimental values. The range- velocity 
and stopping power equations are of great practical inportance 
in a number of fields/ such as in heavy ion acceleration, in 
developing fission electric cells for power generation and in 
the field of radiation bilogy, since these equations allow quan- 
titative prediction of the amoxmt of energy deposited by fast 
moving ions and fission fragments in any conpound media of known 
molecular corrposition like proteins. Heavy ion beams are poten- 
tially the most suitable form of "radiation" for delivering a 
highly localized dose to any deep seated region of choice in the 
body. The following are the stopping power and range-velocity 
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equations: 


Completely stripped heavy ions 

(A) V (2Z^V 
X In 


(a) V)>fz^2^^c 


= 63.65 — ■ log ( 


dx 


V 


11.39 V 

2VV I 
1 o 


A 

[R 2 -R 1 ] = (3.494x10-^) ^ 


A jV3 


Ei{ln C^) -Ei(lnC^)J 


where Cj, or in general C is given by 


C = 


11.39 j 2/3 

. ^V ^i , J 


(b) ifz^^v} 


,1/4 3/4 „ 

'1 ^2 o 


dE r 

~ = 13.79 ^ 3(S„-2) + -yr 

dx ^ ^j2 \ 2 ■ Z 

^'2 


i~) + 2.303 




2Vf (Z^) 


V 


'^6 log 2 V + SCZ^-Z) log y (2 -2) ”^2 2Z^ ^ 

2 o o 2 


(•c) v4' 


1/4 ^ 3/4 ^ 


Zi- - Z^ 


f2 


f (Z ) Z ' 
= 12.68 a ' " ' ^ ~" T" 

dx Ao V 


2 f h-t] 


V ,“i 


V .-1/3 ( , 

„ Of , ry 


+ {2 Z 


1 V 


) 
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(&.) vy 2Z V 


1 o 


d? ^ g3^g5 ^!_|2 ( 1 U ?9 V^ ^ 


dx 




r \ ^ ‘‘~^l 

i R_-Rt [=(1.435x10 n — ^ 

r7 2 


1 


7, 7 ^ 


Ei(ln d 2 ) -Si(ln d^) 


where d^/ <^2 general d is given by 

2 


d = V^) 

X 


and holding for 


v> ( 


y 11.39 


1/2 


Partially stripped ions and fission fracnnents 


m 

dx 


f (Zj 


= 1,327 


.5/3 


4.7622k(Z.)/'' + f(Sj V 

I ■'• J 


R = 


A^A^CVVo) 


127.3 f(Z2)|4.7622|f(Z^)<^^^ + fCZ^)' 

. . dS . . .. ^ MeV cm^ , 

In all these equations — is in units of , while 

^ dx mg 

2 

(R^-R^) or R is in units of mg/cm . A^, A^ and Z^, are the 

mass numbers and the atomic nuxribers of the ion and the medium 


respectively. is the initial velocity of the ion, V is the 

velocity at which the stopping power is being calculated, and 

V = e^/la(2.18xlO® cn/sec) . V., V and are to be expressed 

G i O 

R — • 

in units of 10 crr/sec. I is the mean ionization potential of 
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the medium in eV and fiz^) or f iz^) or in general f{z) is 
given by 


f(Z) = 0.28 Z^^^ 

for Z j^45 .5 

f(Z) = z^/^ 

for Z.^45.5 


The equations given above are applicable in the case of compound 
media if A. and f(Z ) are replaced by E.X. A. and 2.X.f(Z.) 

^ w *3. 3L 3. . 3L 3>i. 

til 

respectively, where the conpound is made up or X^. atoms of the i 
atomic species having mass number and atomic number A^ and Z^, 
respectively! 
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